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ABSTRACT 

We report measurements of the cosmic microwave background (CMB) power spectrum from the 
complete 2008 South Pole Telescope (SPT) data set. We analyze twice as much data as the first 
SPT power spectrum analysis, using an improved cosmological parameter estimator which fits multi- 
frequency models to the SPT 150 and 220 GHz bandpowers. We find an excellent fit to the measured 
bandpowers with a model that includes lensed primary CMB anisotropy, secondary thermal (tSZ) and 
kinetic (kSZ) Sunyaev-Zel'dovich anisotropics, unclustered synchrotron point sources, and clustered 
dusty point sources. In addition to measuring the power spectrum of dusty galaxies at high signal-to- 
noise, the data primarily constrain a linear combination of the kSZ and tSZ anisotropy contributions 
at 150 GHz and £ = 3000: 0^^^ + Q.bD^^^^ = 4.5 ± 1.0 ^iK^. The 95% confidence upper limits on 
secondary anisotropy power are D\qqq < 5.3 ^K^ and -Dfooo < 6-5 fiK^. We also consider the potential 
correlation of dusty and tSZ sources, and find it incapable of relaxing the tSZ upper limit. These 
results increase the significance of the lower than expected tSZ amplitude previously determined from 
SPT power spectrum measurements. We find that models including non-thermal pressure support in 
groups and clusters predict tSZ power in better agreement with the SPT data. Combining the tSZ 
power measurement with primary CMB data halves the statistical uncertainty on as- However, the 
preferred value of erg varies significantly between tSZ models. Improved constraints on cosmological 
parameters from tSZ power spectrum measurements require continued progress in the modeling of the 
tSZ power. 

Subject headings: cosmology - cosmology:cosmic microwave background - cosmology: observations - 
large-scale structure of universe 
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1. INTRODUCTION 



Measurements of temperature anisotropy in the cos- 
mic microwave background (CMB) have proven to be 
some of the most pow erful an d robust tests of cosmo- 
logical theory ( Komat su et al.p OlO; Larson et al. |2010 
Reichardt et al.||2009i lUunkley et al.,i20iul. i'he large- 



and intermediate-scale anisotropy is dominated by the 
signal from the primary anisotropy of the CMB. On 
smaller scales, the primary anisotropy is exponentially 
suppressed by Silk damping (Silk 1968); fluctuations on 
these scales are dominated by foregrounds and secondary 
anisotropics. These secondary anisotropics are produced 
by interactions between CMB photons and matter that 
lies between the surface of last scattering and the ob- 
server. Only recently have experiments reached sufficient 
sensitivity and resolution to explore the cosmological in- 
formation encoded in the secondary CMB anisotropies. 

The largest of these secondary anisotropies is expected 
to be the Suny aev-Zel'dovich (SZ) effect (Sunyaev & 
Zerdovich||1972 ) , which consists of two components: the 
thermal SZ (tSZ) and kinetic SZ (kSZ) effects. The tSZ 
effect occurs when CMB photons are inverse Compton 
scattered by hot electrons in the gravitational poten- 
tial well of massive galaxy clusters. The resulting spec- 
tral distortion of the CMB blackbody spectrum creates 
a decrement in intensity at low frequencies, an incre- 
ment at high frequencies, and a null near 220 GHz. At 
150 GHz, the anisotropy induced by the tSZ effect is ex- 
pected to dominate over the primary CMB fluctuations 
at multipoles I > 3500. 

The anisotropy power due to the tSZ effect depends 
sensitively on the normalization of the matter power 
spectrum, as parametrized by the RMS of the mass dis- 
tribution on 8h~^ Mpc scales, erg. Measurements of the 
tSZ power spectrum have the potential to provide inde- 
pendent constraints on cosmological parameters such as 
cTg and improve the precision with which they are de- 
termined. In principle, the tSZ power spectrum from 
large surveys can also be used to constrain non-standard 
cosmological models, for example, placing l imits on the 
range of allowed early dark energy models (Alam et al. 
20l0| . 

The tSZ power spectrum is challenging to model accu- 
rately because it includes significant contributions from 
galaxy clusters that span a wide range of mass and 
redshift. Modeling uncertainties arise from both non- 
gravitational heating effects in low-mass clusters and the 
limited body of observational data on high-redshift clus- 
ters. As a result, it has been difficult to accurately model 
the shape and amplitude of the tSZ power spectrum ex- 
pected for a given cosmology. Recent models that vary in 
their treatmen t of cluster gas physics differ in am plitude 
by up to 50% ( |Shaw et al j|2010l ' |Trac eraL][2010l) 

The kSZ effect occurs wlien CMB photons are JJoppler 
shifted by the bulk velocity of electrons in intervening 
gas. The kSZ power spectrum is not weighted by the gas 
temperature and is, therefore, less sensitive to the non- 
linear effects that complicate tSZ models. It does, how- 
ever, depend on the details of reionization, which are not 
yet well understood. In the standard picture of reion- 
ization, ionized bubbles form around the first stars or 
quasars, eventually merging and leading to a fully ion- 
ized Universe. The ionized bubbles will impart a Doppler 



shift on scattered CMB photons proportional to the bub- 
ble velocity. This so-called "patchy" reionization kSZ 
signal depends on t he duration of re ionizati on as well 



as the bubble sizes ( Gruzinov fc Hu, 1998; Knox et al. 
19981. A detection oi the kSZ power spectrum or upper 



limit on its amplitude can, in principle, lead to interest 



ing co nstraints on the epoch of reionization ( Zahn et al. 



20051 



The firs t SPT power spectr um results were recently re- 
ported by Lueker et al. (2010 hereafter LIO). This power 
spectrum was produced with 150 and 220 GHz data 
taken by the South Pole Telescope (SPT) on 100 deg^ 
of sky. The two frequencies were combined to remove 
foregrounds with a dust-like spectrum, resulting in the 
detection of a linear combination of kSZ and tSZ power. 
LIO reported the best-fit amplitude of the SZ power 
{Dt = C 1,1(1 l)/27r) at a multipole of ^ = 3000 to be 
^3000 + 0-46At^o^ = 4.2 ± 1.5AiK2. This amplitude was 
less than expected for many models, implying that either 
CTg was at the low end of the range allowed by other mea- 
surements, or the models were over-predicting the tSZ 
power . The SPT 150 GHz bandpowers were used by Hall| 
et al. (2010 hereafter HIO) to set a 95% CL upper limit 
of 13 /xK^n the kSZ power at £ = 3000. 

The Atacama Cosmology Telescope (ACT) collabora- 
tion has also measured the liigh-^ power spec trum (|Das 



et al. 



as also measured tne liigli-^ power spec trum ([JJas 
2010[ ). [Punkley et al.| ( |20"To] ) use the |Das etal 



( [2010 1 bandpowers to measure the sum of the tSZ and 



kSZ power at £ = 3000 and 148 GHz to be 6.8 ± 2.9 /.tK^. 
The SPT and ACT bandpowers and resulting constraints 
on SZ power are consistent within the reported uncer- 
tainties. 

Since the publication of LIO, several new models for 



the tSZ power spectrum have been published ( Trac et al. 
2010| [Shaw et al.||2010[ [Battagha et aL||2010[ ) . Most tea- 



ture similar angular scale dependencies (indistinguish 
able with current data); however, the predicted ampli 
tude of the SZ signal varies considerably between models. 
The amplitude is generally reduced as the amount of star- 
formation, feedback, and non-thermal pressure support 
in clusters increases. These effects can re duce the pre- 
dicted tSZ power by 50% compared to thelSehgal et al 



( 2010[ SIO) model used by LIO, and greatly reduce the 



tension between the measured and predicted tSZ power. 

In addition to the tSZ and kSZ signal, the measured 
power at these frequencies and angular scales includes 
contributions from several significant foregrounds. The 
power from dusty star-forming galaxies (DSFGs) has 
both a Poisson and clustered component, with distinct 
angular scale dependencies. Measurements of DSFG 
power at multiple millimeter wavelengths can be com- 
bined to constrain the DSFG spectral index. HIO con- 
strained the spectral index of the Poisson term to be 
OLy = 3.86 ± 0.23 and the index of the clustered term to 
be Oic = 3.8 ± 1.3 using fits to the LIO single-frequency 
bandpowers. HIO also interpret the implications of the 
spectral index constraint on the dust te mperature and 
redshi ft distribution of the dusty galaxies. [Punkley et al.| 
(20101 also detect significant power attributed to clus- 



tereT DSFGs in the ACT data. They find a preferred 
spectral index of 3.69 ± 0.14, consistent with the HIO 
result. 

In this work, we make two key improvements upon the 
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first SPT power spectruni release (LIO and HIO). First, 
we include both fields observed by SPT in 2008, doubling 
the sky coverage and reducing the bandpower uncertain- 
ties by 30%. Second, we test cosmological models with 
a true multi-frequency analysis of the bandpowers, prop- 
erly accounting for the multi-frequency covariance matrix 
and frequency dependence of each component to esti- 
mate cosmological parameters. Including the additional 
frequency information in the parameter estimation leads 
to improved model constraints. We present the minimal 
model required to explain the data, and then explore a 
number of extensions to this minimal model such as the 
possibility of the tSZ signal being correlated with dust 
emission. None of these extensions significantly improve 
the fits or change the general conclusions concerning the 
amplitudes of the kSZ and tSZ power. 

We describe the instrument, observations, beams, and 
calibration strategy in ^|2j The time-ordered data filter- 
ing and map-making algorithm is outlined in ^ along 
with the procedure to derive bandpowers from maps. 
The results of tests for systematic errors applied to the 
SPT data are discussed in f|4j The bandpowers are given 
in S}5j and the model is presented in Sj6] The fitted pa- 
rameters and their cosmological interpretation are given 
in fjTland |8j respectively. We summarize our conclusions 
in |9] 

2. INSTRUMENT AND OBSERVATIONS 

The SPT is an off-axis Gregorian telescope with a 
10 meter diameter primary mirror located at the South 
Pole. The receiver is equipped with 960 horn-coupled spi- 
derweb bolometers with superconducting transition edge 
sensors. In 2008, the focal plane included detectors at 
two frequency bands centered at approximately 150 and 
220 GHz. The telescope and rec ei ver are discussed in 
mor e detail in |Ruhl et "aD ( |2004| ), |Padin et al^ ( |2008[ ), 
and Carlstrom et ai. 1 1UT^ . 

In this work, we use data at 150 and 220 GHz from two 
--100 deg2 fields observed by SPT in the 2008 austral 
winter. The fields are centered at right ascension (RA) 
5^30™, dechnation (decl.) -55° (J2000) (henceforth the 
5^30™ field) and RA 23'^30'", decl. -55° (the 23*^30™ 
field) . The locations of the fields were chosen for overlap 
with the optical Blanco Cosmology Survey (BC Srjand 
low dust em ission observed by IRAS at 100 fim (IScmegel 
et al. 1998 1. After data quality cuts, a total of 1100 hours 
of observations are used in this analysis. The final map 
noise is 18 /iK-arcmirj^ at 150 GHz and 40 /iK-arcmin at 
220 GHz. These data include the majority of the sky 
area observed in 2008. In 2009 and 2010, the SPT has 
been used to survey an additional 1300 deg^ to the same 
depth at 150 GHz. 

The fields were observed with two different scan 
strategies. The 5^30"^ scan strategy involves constant- 
elevation scans across the 10° wide field. After each scan 
back and forth in azimuth across the field, the telescope 
executes a 0.125° step in elevation. A complete set of 

http://cosmology.illinois.edu/BCS 

Throughout this work, the unit K refers to equivalent fluctu- 
ations in the CMB temperature, i.e., the temperature fluctuation 
of a 2.73 K blackbody that would be required to produce the same 
power fluctuation. The conversion factor is given by the derivative 
of the blackbody spectrum, evaluated at 2.73 K. 



scans covering the entire field takes approximately two 
hours, and we refer to each complete set as an observa- 
tion. 

The 23*^30™ field was observed using a lead-trail scan 
strategy. Two consecutive observations each map ap- 
proximately half the field. Due to the Earth's rotation, 
both halves of the field are observed at the same range of 
azimuth angle. This would allow the removal of ground- 
synchronous signals if detected; however, we see no evi- 
dence for such a signal at the angular scales of interest. 
A complete observation of one half-field takes approx- 
imately 40 minutes. The two halves are generally ob- 
served directly after one another, though for this analysis 
we also include a small number of lead trail pairs (31 out 
of 480) which were not observed sequentially, due to poor 
weather or other interruptions. The requirement that 
both halves of the 23^30™ field are observed at the same 
azimuth leads to a larger elevation step (0.268°) and 
therefore less uniform coverage across the map. We ame- 
liorate this non-uniform coverage in the cross-spectrum 
analysis by combining two lead-trail map pairs into one 
'observation' unit. The lead-trail pairs are chosen to to 
maximize coverage uniformity in the resulting map and 
minimize temporal offsets. 

2.1. Beam Functions 

The power spectrum analysis presented here depends 
on an accurate measurement of the beam function, which 
is the azimuthally averaged Fourier transform of the 
beam map. Due to the limited dynamic range of the 
detectors, the SPT beams for the 2008 observing sea- 
son were measured by combining maps of three sources: 
Jupiter, Venus, and the brightest point source in each 
field. Maps of Venus are used to stitch together the outer 
and inner beam maps from Jupiter and the point source, 
respectively. Maps of Jupiter at radii > 4' are used to 
constrain a diffuse, low-level sidelobe that accounts for 
roughly 15% of the total beam solid angle. The beam 
within a radius of 4' is measured on the brightest point 
source in the field. The effective beam is slightly en- 
larged by the effect of random errors in the pointing re- 
construction. We include this effect by measuring the 
main beam from a source in the final map. The main- 
lobe beam is approximately fit by 2D Gaussians with 
full width at half-maximum (FWHM) equal to 1.15' at 
150 GHz and 1.05' at 220 GHz. We have verified that the 
measured beam is independent of the point source used, 
although the signal-to-noise drops for the other (dim- 
mer) sources in either field. The beam measurement is 
similar to that described in LIO, although the radii over 
which each source contributes to the beam map have 
changed. However, the estimation of beam uncertainties 
has changed. We now determine the top three eigenvec- 
tors of the covariance matrix for the beam function at 
each frequency, and include these as beam errors in the 
parameter fitting. These three modes account for > 99% 
of the beam covariance power. The beam functions and 
the quadrature sum of the three beam uncertainty pa- 
rameters are shown in Figure [T] 

2.2. Calibration 

The absolute calibration of the SPT data is based 
on a comparison of the CMB power at degree scales in 
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Fig. 1. — Left axis: The measured SPT beam functions at 
150 GHz (black line) and 220 GHz (blue line). Right axis: The 
fractional beam uncertainties at 150 GHz (black dashed line) 
and 220 GHz (blue dashed line). The beam uncertainty is 
parametrized by a three component model, the quadrature sum 
of which is plotted here. 

the WMAP 5-year maps with dedicated SPT cahbration 
maps. These cahbration maps cover four large fields to- 
taling 1250 deg^ which were observed to shallow depth 
during the 2008 season. Details of the cross-calibration 
with WMAP are given in LIO, and remain essentially 
unchanged in this work. Although the beam functions 
changed slightly at £ > 1000 as a resu lt of the improved 
beam measurement described in §2.11 this change had a 
negligible effect at the angular scales iased in the WMAP 
calibration. The calibration factors are therefore un- 
changed from LIO; however as discussed above, we have 
slightly changed the treatment of beam and calibration 
errors, and are no longer folding part of the beam uncer- 
tainty into a calibration uncertainty. 

We estimate the uncertainty of the temperature cali- 
bration factor to be 3.5%, which is slightly smaller than 
the value of 3.6% presented in LIO. We applied this 
calibration to the 220 GHz band by comparing 150 GHz 
and 220 GHz estimates of CMB anisotropy in the survey 
regions. This internal cross-calibration for 220 GHz is 
more precise than, but consistent with, a direct absolute 
calibration using observations of RCW38. We estimate 
the 220 GHz temperature calibration uncertainty to be 
7.1%. Because the 220 GHz calibration is derived from 
the 150 GHz calibration to WMAP, the calibration uncer- 
tainties in the two bands are correlated with a correlation 
coefficient of appro xima tely 0.5. 

As discussed in |7.1[ our MCMC chains also include 
two parameters to adjust the overall temperature cal- 
ibration in each SPT band. The resulting calibration 
factors are listed in the Table [T] comments. 

3. ANALYSIS 

In this section, we present an overview of the pipeline 
used to process the time-ordered data (TOD) to band- 
powers. The method closely follows the approach used by 
LIO, and we refer the reader to LIO for complete details 
of the method. We highlight any differences between this 
analysis and that work. 

Given the small field sizes, we use the flat-sky approx- 
imation. We generate maps using the oblique Lambert 
equal-area azimuthal projection and analyze these maps 
using Fourier transforms. Thus, the discussion of filtering 
and data-processing techniques refer to particular modes 
by their corresponding angular wavenumber k in radians 
with \k\=l 



3.1. TODs to Maps 

Each detector in the focal plane measures the CMB 
brightness temperature plus noise, and records this mea- 
surement as the TOD. Noise in the TOD contains con- 
tributions from the instrument and atmosphere. The 
instrumental noise is largely uncorrelated between detec- 
tors. Conversely, the atmospheric contribution is highly 
correlated across the focal plane, because the beams of 
individual detectors overlap significantly as they pass 
through turbulent layers in the atmosphere. 

We bandpass filter the TOD to remove noise outside 
the signal band. The TOD are low-pass filtered at 12.5 
Hz {£ ^ 18000) to remove noise above the Nyquist fre- 
quency of the map pixelation. The TOD are effectively 
high-pass filtered by the removal of a Legendre polyno- 
mial from the TOD of each detector on a by-scan basis. 
The order of the polynomial is chosen to have the same 
number of degrees of freedom (dof) per unit angular dis- 
tance (^1.7 dof per degree). The polynomial fit removes 
low frequency noise contributions from the instrument 
and atmosphere. 

Correlated atmospheric signals remain in the TOD af- 
ter the bandpass filtering. We remove the correlated sig- 
nal by subtracting the mean signal over each bolometer 
wedge ^ at each time sample. This filtering scheme is 
slightly different from that used by LIO, although nearly 
identical in effect at 150 GHz. LIO subtracted a plane 
from all wedges at a given observing frequency at each 
time sample. Since SPT had two wedges of 220 GHz 
bolometers and three wedges of 150 GHz bolometers in 
2008, the LIO scheme filters the 220 GHz data more 
strongly than the 150 GHz data. The wedge-based com- 
mon mode removal implemented in this work results 
in more consistent filtering between the two observing 
bands. 

The data from each bolometer is inverse noise weighted 
based on the calibrated, pre-filtering detector power 
spectral density. We bin the data into map pixels based 
on pointing information, and in the case of the 23^30™ 
field, combine two pairs of lead and trail maps to form 
individual observations as discussed in [|2j A total of 
300 and 240 individual observation units are used in the 
subsequent analysis of the 5*^30™ and 23*^30™ fields, re- 
spectively. 

3.2. Maps to Bandpowers 

We use a pseudo-Q method to estimate the bandpow- 
ers. In pseudo-Cf methods, bandpowers are estimated 
directly from the Fourier transform of the map after cor- 
recting for effects such as TOD filtering, beams, and fi- 
nite sky coverage. We p rocess the data usin g a cross 



spectrum b ased analysis f Polenta et al.| 2005 



et al. 



Tristram 
Beam and 



2005[ ) in order to eliminate noise bias 
filtering ettects are corrected for according to the fo rmal- 
ism in the MASTER algorithm ( |Hivon et aL]|2002[ ). We 
report the bandpowers in terms ot 'L>£^ where 



27r 



(1) 



^® The SPT array consists of 6 pie shaped bolometer wedges 
of 160 detectors, of which five wedges are used in this analysis. 
Wedges are configured with a set of filters that determine their 
observing frequency (e.g. 150 or 220 GHz). 
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The first step in the analysis is to calculate the Fourier 
transform irS^'^'^'^ of the map for each frequency vi and 
observation A. All maps of the same field are apodized by 
a single window that is chosen to mask out detected point 
sources and avoid sharp edges at the map borders. After 
windowing, the effective sky area used in this analysis is 
210 deg^ . Cross-spectra are then calculated for each map- 
pair on the same field; a total of 44850 and 28680 pairs 
are used in the S'^SO^and 23'^30™fields, respectively. We 
take a weighted average of the cross-spectra within an 
£-bin, b. 



v,-KUj,AB _ I k(k + 1) , _ {vj,B)* 



2tt 



(2) 



keb 



As in LIO, we find a simple, uniform selection of modes at 
kx > 1200 is close to the optimal mode-weighting. With 
the adopted flat-sky approximation, £ = k. 

The bandpowers above are band-averaged pseudo-C^s 
that can be related to the true astronomical power spec- 
trum D by 



D 



fiXUj ,AB 



(3) 



The K matrix accounts for the effects of timestream fil- 
tering, windowing, and band-averaging. This matrix can 
be expanded as 



Kt^"' ^Pbk{MkAw]F;:r 



'Bk'B'k' ) Qk'b' 



(4) 



Here Qk h and Phk are the binning and re-binning op- 



erators (|Hivon et al. 2002|). _Bj^* is the beam function 



for frequency v^, and Fj^^ is the fc-dependent transfer 
function which accounts for the filtering and map-making 
procedure. The mode-mixing due to observing a finite 
portion of the sky is represented by M/cfc/[W], which is 
calculated analytically from the known window W. 

3.2.1. Transfer Function 

We calculate a transfer function for each field and ob- 
serving frequency. The transfer functions are calculated 
from simulated observations of 300 sky realizations that 
have been smoothed by the appropriate beam. Each sky 
realization is a Gaussian realization of the best-fit lensed 
WMAP7 ACDM model plus a Poisson point source con- 
tribution. The Poisson contribution is selected to be 
consistent with the results of HIO, and includes radio 
source and DSFG populatio ns. The amplitude of the ra- 
dio source term is set by the de Zotti et al. ( 2005 1 model 
source counts at 150 GHz with an assumed spectral in- 
dex of Ur = —0.5. The DSFG term has an assumed spec- 
tral index of 3.8. The effective point source powers are 
^150x150 ^ 7 5 X 10-6 ^K2, 0''22o ^ 2.3 X 10-5 /xK2, 

and C'|20'^220 ^ 7 g ^ lO-^/xK^. The sky reahzations 
are sampled using the SPT pointing information, filtered 
identically to the real data, and processed into maps. 
The power spectrum of the simulated maps is compared 
to the known input spectr um, (7'^.th°°i'y to ca lculate the 
effective transfer function ( Hivon et al. 2002 1 in an iter- 
ative scheme, 



^k — ^k 



MC 



Mkk'i'k' ^k' ^k' 



(5) 



Here 102 = J cfcW^ is a normalization factor for the area 
of the window. We find that the transfer function esti- 
mate converges after the first iteration and use the fifth 
iteration. For both fields and both bands, the result is a 
slowly varying function of £ with values that range from 
a minimum of 0.6 at ^ = 2000 to a broad peak near 
£ = 5000 of 0.8. 

3.2.2. Bandpower Covariance Matrix 

The bandpower covariance matrix includes two terms: 
variance of the signal in the field (cosmic variance), 
and instrumental noise variance. The signal-only Monte 
Carlo bandpowers are used to estimate the cosmic vari- 
ance contribution. The instrumental noise variance is 
calculated from the distribution of the cross-spectrum 



bandpowers D 



lyiXUj ,AB 



between observations A and B, 



as described in LIO. We expect some statistical uncer- 
tainty of the form 



(6) 



in the estimated bandpower covariance matrix. Here, 
riobs is the number of observations that go into the co- 
variance estimate. This uncertainty on the covariance is 
significantly higher than the true covariance for almost 
all off-diagonal terms due to the dependence of the un- 
certainty on the (large) diagonal covariances. We reduce 
the impact of this uncertainty by "conditioning" the co- 
variance matrix. 

How we condition the covariance matrix is determined 
by the form we expect it to assume. The covariance ma- 
trix can be viewed as a set of nine square blocks, with 
the three on-diagonal blocks corresponding to the covari- 
ances of a 150x150, 150x220, or 220 x220 GHz spectrum. 
Since the bandpowers reported in Tables [T] and |4] are ob- 
tained by first computing power spectra and covariance 
matrices for bins of width A£ = 100 with a total of 80 ini- 
tial £-bins, each of these blocks is an 80 x 80 matrix. The 
shape of the correlation matrix in each of these blocks 
is expected to be the same, as it is set by the apodiza- 
tion window. As a first step to conditioning the covari- 
ance matrix, we calculate the correlation matrices for the 
three on-diagonal blocks and average all off-diagonal el- 
ements at a fixed separation from the diagonal in each 
block, 

ki — k2—k — k' p^, \ \ 

^kk' — ^ ■ \' ) 

A^ki —k2 — k—k' 

This averaged correlation matrix is then applied to all 
nine blocks. 

The covariance matrix generally includes an estimate 
of both signal and noise variance. However, the covari- 
ance between the 150 x 150 and 220 x 220 bandpowers 
is a special case: since we neither expect nor observe 
correlated noise between the two frequencies in the sig- 
nal band, we include only the signal variance in the two 
blocks describing their covariance. 

3.2.3. Combining the Fields 

We have two sets of bandpowers and covariances — one 
per field — which must be combined to find the best esti- 
mate of the true power spectrum. We calculate near- 
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optimal weightings for this combination at each fre- 
quency and i-h'm using the diagonal elements of the co- 
variance matrix, 



(8) 

This produces a diagonal weight matrix, w, for which the 
bandpowers can be calculated: 



Db =!Li^c 

c 



(9) 



where C, denotes the field. The covariance will be 

(l/i Xl'j ) X (I'm X I/„) 

^bb' ~ 

E( ViXVi (l/i XI'J ) X (I'm X i/„xi/„\ /inN 

c 

After combining the two fields, we average the bandpow- 
ers and covariance matrix into the final ^-bins. Each 
initial A£ — 100 sub-bin gets equal weight. 

4. JACKKNIFE TESTS 

We apply a set of jackknife tests to the data to search 
for possible systematic errors. In a jackknife test, the 
data set is divided into two halves based on features of 
the data associated with potential sources of systematic 
error. After differencing the two halves to remove any as- 
tronomical signal, the resulting power spectrum is com- 
pared to zero. Significant deviations from zero would 
indicate a systematic problem or possibly a noise mises- 
timate. Jackknives with a cross-spectrum (as opposed to 
an auto-spectrum) estimator are less sensitive to small 
noise misestimates since there is no noise bias term to be 
subtracted. We implement the jackknives in the cross- 
spectrum framework by differencing single pairs of obser- 
vations a nd a pplying the cross-spectrum estimator out- 
lined in ^3.2 to the set of differenced pairs. In total, 
we perform three jackknife tests based on the observing 
parameters: time, scan direction, and azimuthal range. 

The data can be split based on the time of observation 
to search for variability in the calibration, beams, detec- 
tor time constants, or any other potentially time vari- 
able aspect of the observations. The "first half - second 
half" jackknife probes variations on month time scales. 
Results for the "first half - second half" jackknife are 
shown in the top panel of Figure [2j We note that in the 
23*^30™ field, the combination of individual field observa- 
tions into four-map units, discussed in section [2| results 
in a small number of constituent maps (9 out of 960) that 
are grouped into the wrong side of this split. This does 
not affect the 5^30'" field, where the observation unit is 
a single map. 

The data can also be split based on the direction of 
the scan in a "left -right" jackknife (panel 2 of Figure [2]) . 
We would expect to see residual power here if the detec- 
tor transfer function has been improperly de-convolved, 
if the telescope acceleration at turn-arounds induces a 
signal through sky modulation or microphonics, or if the 
wind direction is important. 

Sidelobe pickup could potentially introduce spurious 
signals into the map from features on the ground. In 
deep coadds of unfiltered data, we observe features with 
scales of several degrees that are fixed with azimuth and 



are presumably caused by ground pickup. This ground 
pickup is significantly reduced by the wedge common- 
mode removal, but could still exist at low levels. We 
split the data as a function of azimuth based on the ob- 
served RMS signal on these large scales. This split is 
different than the choice made in LIO, which was based 
upon the distance in azimuth between an observation and 
the closest building. The azimuthal jackknife is shown in 
the third panel of Figure [2] 

We calculate the with respect to zero for each 
jackknife over the range i € [2000, 9500] in bins with 
A£ = 500. The probability to exceed the measured 
(PTE) for the three individual jackknives is 5%, 66%, 
and 15% at 150 GHz and 99%, 21%, and 98% at 220 
GHz for the "first half - second half," "left- right," and 
"azimuth" jackknives, respectively. The combined PTE 
for the three jackknives is 10% for the 150 GHz data, 96% 
for the 220 GHz data and 55% for the combined set of 
both frequencies. We thus find no evidence for system- 
atic contamination of the SPT band powers. 

5. POWER SPECTRA 

The bandpowers presented in Table [l] and plotted in 
Figure [3] are the result of applying the analysis described 
in f|| to the 210 deg^ observed by SPT in 2008. The 
errors shown in Table [T] are the diagonal elements of the 
covariance matrix. The full covariance matrix includin; 
off-diagonal elements can be found at the SPT websit' 
along with the window functions describing the transfor 
mation from a theoretical spectrum to these bandpowers 



(Knox 



1999) 



In these power spectra, we have measured the mm- 
wave anisotropy power at £ > 2000 with the highest sig- 
nificance to date. Anisotropy power is detected at 92 a 
at 150 GHz and 74 cr at 220 GHz. Below £ - 3000, the 
bandpowers refiect the Silk damping tail of the primary 
CMB anisotropy. At ^ > 3000, the bandpowers follow 
the Di (X i'^ form expected for a Poisson point source 
distribution. The measured spectral index (a ~ 3.6 
where Si, cx v") points to a dust-like spectrum for the 
point sources indicating that the DSFGs are the domi- 
nant point source population. However, the CMB and 
Poisson terms do not fully explain the data. Adding two 
free parameters representing the tSZ and effects of point 
source clustering improves the best-fit model likelihood, 
£, by Aln£ — 65. Using a combination of the 150 GHz 
auto-spectrum, 150 x 220 cross-spectrum and 220 GHz 
auto-spectrum, we can separate these foregrounds and 
secondary anisotropics; however, substantial degenera- 
cies remain. Future SPT data analyses will include a 
third frequency band, which we expect will significantly 
reduce these degeneracies. 

6. COSMOLOGICAL MODEL 

We fit the SPT bandpowers to a model including 
lensed primary CMB anisotropy, secondary tSZ and kSZ 
anisotropics, galactic cirrus, an unclustcred population of 
radio sources and a clustered population of DSFGs. In 
total, this model has six parameters describing the pri- 
mary CMB anisotropy, up to two parameters describing 
secondary CMB anisotropics, and up to eleven parame- 
ters describing foregrounds. 

'^'^ http:/ /pole. uchicago.edu/public/data/shirokofT10/ 
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Fig. 2.— Jackknives for the SPT data set at 150 GHz (blue circles) and 220 GHz (black diamonds). For clarity, the 220 GHz 
jackknives have been shifted to the right by A£ = 100. Top panel: Bandpowers of the "first half - second half" jackknife compared to the 
expected error bars about zero signal. Disagreement with zero would indicate either a noise misestimate or a time-dependent systematic 
signal. Second panel: Power spectrum of the left-going minus right-going difference map. This test yields strong constraints on the accuracy 
of the detector transfer function deconvolution and on possible directional systematics. Third panel: Bandpowers for the difference map 
when the data are split based on the observed very large scale {£ < 100) ground pickup as a function of azimuth. Signals fixed in azimuth, 
such as ground pickup on smaller scales, would produce non-zero power. We see no evidence for ground-based pickup across this £-range. 
The cumulative probability to exceed the observed in these three tests at 150 and 220 GHz is 55%. Bottom panel: The un-differenced 
SPT power spectra at each frequency for comparison. 

sons to expect contributions to the power spectra from 
each of the components outhned above, a good fit to 
the current SPT data can be obtained with a much 
simpler, four-parameter extension to the primary CMB 
anisotropy. This minimal model contains four free pa- 
rameters: the tSZ amplitude, the Poisson and clustered 
DSFG amplitudes, and a single spectral index which sets 
the frequency dependence of both DSFG terms. All other 
parameters are tightly constrained by external data. A 
similar model was used in recent ACT analysis ( Dunkley 



Parameter constraints are calc ulated using the publicly 
available CosmoMCj^ package ([Lewis fc Bridle] |2002| . 
We have added two modules to nandle the high-* data; 
one to model the foregrounds and secondary anisotropics 
and one to calculate the SPT likelihood function. These 
extensions, with some differences in the tre atment of the 
second ary anisotropy, are also discussed in |Millea et al.| 
(2010). These modules and instructions for compiling 



them are available at the SPT website 1^ 
Although there are observational and theoretical rea- 

http://cosmologist.info/cosmomc 
■^^ http: / /polo. uchicago.edu/public/data/shirokofflO/ 



et al.| |2010). We adopt this four-parameter extension as 
the baseline model for this work. 
In this baseline model, several simplifications are made. 
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TABLE 1 
Single-frequency bandpowers 









150 GHz 


150 X 220 GHz 


220 GHz 


I range 


^Gir 


b (/iK^) 


a (mK^) 


b {^,K^) 


o (mK^) 


b (mK2) 


o (mK^) 


2001 - 


2200 


2056 


242.1 


6.7 


248.8 


8.3 


295.9 


14.9 


2201 - 


2400 


2273 


143.2 


4.2 


154.7 


5.6 


201.5 


11.7 


2401 - 


2600 


2471 


109.3 


3.2 


122.1 


4.5 


193.5 


11.0 


2601 - 


2800 


2673 


75.9 


2.6 


102.8 


4.1 


172.1 


10.4 


2801 - 


3000 


2892 


60.2 


2.3 


80.4 


3.7 


179.3 


11.1 


3001 - 


3400 


3184 


47.5 


1.5 


73.5 


2.5 


169.7 


8.2 


3401 - 


3800 


3580 


36.9 


1.6 


69.5 


2.7 


180.7 


9.2 


3801 - 


4200 


3993 


36.7 


1.8 


81.0 


3.3 


240.1 


11.6 


4201 - 


4600 


4401 


38.5 


2.2 


81.5 


3.9 


226.6 


12.6 


4601 - 


5000 


4789 


39.3 


2.7 


96.9 


4.7 


276.9 


15.2 


5001 - 


5900 


5448 


49.2 


2.5 


122.3 


4.2 


361.3 


13.3 


5901 - 


6800 


6359 


60.8 


3.9 


158.7 


6.2 


457.4 


20.0 


6801 - 


7700 


7256 


89.1 


6.0 


173.9 


9.5 


488.8 


28.4 


7701 - 


8600 


8159 


81.7 


9.5 


229.2 


14.2 


653.0 


42.3 


8601 - 


9500 


9061 


131.9 


14.4 


309.1 


21.9 


895.2 


66.0 



Note. — Band multipole range and weighted value laS, bandpower b, and uncertainty a for the 150 GHz auto-spectrum, 
cross-spectrum, and 220 GHz auto-spectrum of the SPT fields. The quoted uncertainties include instrumental noise and the 
Gaussian sample variance of the primary CMB and the point source foregrounds. The sample variance of the SZ effect, beam 
uncertainty, and calibration uncertainty are not included. To include the preferred calibration from the MCMC chains, these 
bandpowers should be multiplied by 0.92, 0.95, and 0.98 at 150 GHz, 150 x 220 GHz, and 220 GHz respectively (See j |7T| ) Beam 
uncertainties are shown in Figure n] and calibration uncertainties are quoted in j ]2.2[ Point sources above 6.4mJy at 150 GHz 
have been masked out in this analysis. This flux cut substantially reduces the contribution of radio sources to the bandpowers, 
although DSFGs below this threshold contribute significantly to the bandpowers. 



The distributions of spectral indices are set to delta func- 
tions {fjx = 0). The radio source spectral index is fixed 
to Ur — —0.53. The amplitudes of the kSZ, radio source, 
and cirrus power spectra are fixed to expected values. 
The tSZ-DSFG correlation is set to zero. We discuss 
the priors placed on these components in the following 
subsections. 

To compare the SPT data to these cosmological mod- 
els, we combine the SPT high-£ bandpowers with low-£ 



(iLarson et al.||2010; Reichardt et al. 


,2009; Brown et al. 


2UUyp. Uurmg the preparation ot this paper, the ACT 



it was not until after the model fitting runs were com- 
pleted that the window functions were made available. 
For this reason, we do not include the ACT bandpow- 
ers in the fits although they are the best current con- 
straints for I e [1600,2000]. We expect that including 
the ACT data would have little impact on the results. 
The \ow-l CMB bandpowers tightly constrain the ACDM 
model space. We restrict the ACBAR and QUaD band- 
powers to ^ < 2100 where the details of the secondary 
anisotropics and point source spectra are negligible com- 
pared to the primary CMB anisotropy. These additional 
components are challenging to model between experi- 
ments since they are frequency-dependent and, in the 
case of point sources, mask-dependent. In practice, re- 
moving this restriction does not affect the results since 
the SPT data dominates the constraints above I — 2000. 

6.1. Primary CMB Anisotropy 

We use the standard, six-parameter, spatially flat, 
lensed ACDM cosmological model to predict the primary 



CMB temperature anisotropy. The six parameters are 
the baryon density flbh^, the density of cold dark matter 
Qch^, the optical depth to recombination r, the angu- 
lar scale of the peaks 8, the amplitude of the primordial 
density fluctuations ln[10^°As], and the scalar spectral 
index n^. 

Gravitational lensing of CMB anisotropy by large-scale 
structure tends to increase the power at small angular 
scales, with the potential to influence the high-£ model 
fits. The calculation of lensed CMB spectra out to 
£ — 10000 is prohibitively expensive in computational 
time. Instead, we approximate the impact of lensing by 
adding a fixed lensing template calculated for the best- 
fit WMAP7 cosmology to unlensed CAM^^ spectra. 
The error due to this lensing approximation is less than 
0.5 /iK^ for the allowed parameter space at ^ > 3000. 
This is a factor of 20 (260) times smaller than the sec- 
ondary and foreground power at 150 (220) GHz and 
should have an insignificant impact on the likelihood cal- 
culation. We tested this assumption by including the 
full lensing treatment in a single chain with the base- 
line model described above. As anticipated, we found no 
change in the resulting parameter constraints. 

We find the data from SPT, WMAP7, ACBAR, and 
QUaD considered here prefer gravitational lensing at 
Aln£ = —5.9 (3.4 cr). Most of this detection significance 
is the result of the low-^ data and not the high-^ SPT 
data presented here; without including SPT data in fits, 
lensing is preferred at 3.0 ct. In estimating the preference 
for gravitational lensing, we used the true lensing cal- 
culation instead of the lensing approximation described 

http://camb.info 
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Fig. 3.— Top panel: From left to right, the SPT 150 GHz, 150 X 220 GHz, and 220 GHz bandpowers. Overplotted is the best-fit model 
(red line) with components shown individually. The lensed primary CMB anisotropy is marked by an orange line. The best-fit tSZ 
(purple line) and predicted kSZ (purple dashed line) power spectra are also shown. The predicted radio source term is represented by 
the blue dots. The DSFG Poisson term at each frequency is denoted by the green dashed line and the clustered DSFG component by 
the green dot-dash line. The damping tail of the primary CMB anisotropy is apparent below t = 3000. Above £ = 3000, there is a clear 
excess with an angular scale dependence consistent with point sources. These sources have low fiux (sources detected at >5o" at 150 GHz 
have been masked) and a rising frequency spectrum, consistent with expectations for DSFGs. Bottom panel: Plot of the residual between 
the measured bandpowers and best-fit spectrum. 



above. 

6.2. Secondary CMB Anisotropy 

We parametrize the secondary CMB anisotropy with 
two terms describing the amphtude of the tSZ and kSZ 
power spectra: D\qqq and D^q^q. The model for the SZ 
terms can be written as 



dSZ _ ritsz fvifu2 

^t,vi,U2 ~ -^3000 



^tsz 



■I Vn 



^tSZ 
^3000 



j_ rikSZ 

+ ^3000 



$5 



kSZ 



The assumed cosmological parameters are {fib, ^m, ^A, 
h, Us, cTs) = (0.044, 0.264, 0.736, 0.71, 0.96, 0.80). At 
i = 3000, this model predicts D^f^^^ = 7.5 ^lK^ of tSZ 
power in the SPT 150 GHz band. We use this model in 
all chains where another model is not explicitly specified. 

We also consider tSZ power spectrum mod els reported 
by iTrac et al] |2010|), [B attagUa et~al] ( |2010[ ), and [Shaw 
et al.l (|2U10p. 'I'rac et ai. (2010') followed a procedure 



rf.kSZ ■ 
^3000 



(11) 



Similar to that 



I'rac e 

of yiO, exploring the thermal and kinetic 



SZ power s pectra produ c ed for different input param 



Here, denotes the theoretical model template for 
component X at frequency vq. The frequency depen- 
dence of the tSZ effect is encoded in at the base 
frequency i/q, -D30W),,.o,i'o = ^3000- The kSZ effect has 
the same spectrum as the primary CMB anisotropy, so 
its amplitude is independent of frequency. In this work, 
we set vq to be the effective frequency of the SPT 150 
GHz band (see jj6l|. 



6.2.1. tSZ Power Spectrum 

We adopt four different models for the tSZ power spec- 
trum. Following LIO, we use the power spectrum pre- 
dicted by SIO as the baseline model. SIO combined the 
se mi-analytic model for the intra-cluster medium (ICM) 
of Bode et ah] |2009 ) with a cosmological N-body simula- 
tion to produce simulated thermal and kinetic SZ maps 
from which the template power spectra were measured. 



eters of the Bode et al. (2009) gas model. We adopt 
the nontherrnal20 model (hereaifter the Trac model) pre- 
sented in that work, which differs from the SIO simula- 
tions by having increased star formation, lower energy 
feedback, and the inclusion of 20% non-thermal pres- 
sure support. It predicts a significantly smaller value 
of Df^Q = 4.5 /iK2 when scaled to the SPT 150 GHz 
band. T he second template we consider is that pro- 
duced by Battaglia et al. ([2010) from their Smoothed- 
Particle-Hydrodynamics simulations including radiative 
cooling, star formation and AGN feedback (hereafter the 
Battaglia model). This model predicts -D3000 = 5. 6 /.tK^, 
intermediate between the baseline model and the Trac 
model, a nd peaks at sl ightly higher t than either of those 



models. Shaw et al. ( [2010) investigate the impact of 
cluster astrophysics on the tSZ power spectrum using 
halo model calculations in combination with an analytic 
model for the ICM. We use the baseline model from 
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Fig. 4. — The SPT 150GHz bandpowers (black circles), WMAP7 bandpowers (purple squares), ACBAR bandpowers (green trian- 
gles), QUaD bandpowers (cyan diamonds), and ACT 150GHz bandpowers (orange circles) plotted against the best-fit lensed ACDM 
CMB spectrum. The damping tail of the primary CMB anisotropy is apparent below I = 3000. Above t = 3000, there is a clear excess due 
to secondary anisotropics and residual point sources that has now been measured by both SPT and ACT. Note that the source masking 
threshold in the SPT data (6.4 mjy) is lower than that in the ACT data, so we expect less radio source power at high I. We have multiplied 
the SPT bandpowers by the best-fit calibration of 0.92 as determined in parameter fits. 

which predicts 
The model of 
is also used to rescale all the model 
templates as a tunction of cosmological parameters, as 
described in SjS] 

All four tSZ models exhibit a similar angular scale de- 
pendence over the range of multipoles to which SPT is 
sensitive (see Figure We allow the normalization of 
each model to vary in all chains, and detect similar tSZ 
power in all cases (see Table [6]) . However as we discuss 
in ^ the difference between models is critical in inter- 
preting the detected tSZ power as a constraint on cluster 
physics or erg. 



that work (hereafter the Shaw model) 
Dtsz^ = 4.7 uK^ in the 150 GHz band 
Shaw et al.| ( [20To| 



6.2.2. kSZ Power Spectrum 

We use the SIO (homogeneous reionization) kSZ tem- 
plate as the baseline kSZ model. At £ = 3000, this model 
predicts Dfg^of = 2.05 /iK^ of kSZ power. The kSZ am- 
plitude depends on the details of reionization, and the 
scaling of the kSZ power with cosmological parameters, 
particularly erg, is much weaker than the scaling of the 
tSZ power. We therefore choose to fix the amplitude of 
the kSZ signal to a model value and allow tSZ to vary in- 
depen dently. This treatment differs from [Dunkley et al. 
(2010), which us es a single no rmalization for both SZ 
components, and Millea et al. pOlO) in which a fitting 
function is used to calculate the kSZ power as a function 
of cosmological parameters at each step in the MCMC 
chain. For the current SPT data set, where tSZ and 
kSZ power are largely degenerate, we expect the differ- 
ences in kSZ treatment to be insignificant. This assump- 
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I 

Fig. 5.— Templates used for the tSZ, kSZ, and clustered DSFG 
power discussed in ij6] The top plot shows alternate tSZ templates. 
The black, solid Ime is the (baseline) SIO model. The blue, 
dashed line is the Shaw model. The red, dotted line is the 
Trac model. The teal dot-dash line is the Battaglia model. The 
bottom plot shows both kSZ and clustered DSFG templates. The 
black, solid line is the (baseline) SIO kSZ model. The red, 
dotted line is the patchy kSZ model. The blue, dashed line 
is the (baseline) power-law clustered DSFG template. The teal, 
dot-dash line is the linear-theory clustered DSFG template. The 
clustered DSFG templates have both been normalized to 1 /xK^ at 
I = 3000. 

tion was tested by importance sampling an MCMC chain 
with variable kSZ am plitude according to the scaling at 



3000 described in Millea et al. (2010). As expected 



we find no significant ditterence in ntted parameters 
In ^ we will discuss the impact of two alternate kSZ 
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treatments in addition to the basehne model. First, we 
consider a kSZ template that includes the signal from 
patchy reionization (hereafter the patchy kSZ template). 
This template, which was also used i n LIO, is based upon 
the "FFRT" semi- analytic model of'Za hn et al] ( |2010[ ). 
It was calculated in a 1.5 Gpc/h cosmological column 
where the x- and y-axes correspond to roughly 15 degrees 
on a side and the z-axis corresponds to redshift, with a 
median redshift of 8. The inclusion of patchy reioniza- 
tion increases the kSZ amplitude to at most 3.3 /iK^, and 
changes the shape of D^^^ slightly when compared to 
the fiducial template (see Figure l5| . In the second kSZ 
treatment, we consider the impact of allowing the overall 
kSZ amplitude to vary. In this treatment, we found the 
small difference between the SIO and patchy kSZ tem- 
plate shapes had no impact on the result. We therefore 
use the SIO template for all variable amplitude kSZ con- 
straints. 



6.3. Foregrounds 

At € > 3000, the measured anisotropy is dominated by 
astrophysical foregrounds, particularly bright radio point 
sources. We mask all point sources detected at > 5ct 
(6.4 mjy at 150 GHz) when estimating the power spec- 
trum. In order of importance after masking the brightest 
point sources, we expect contributions from DSFGs (in- 
cluding their clustering), radio sources (unclustered), and 
galactic cirrus. We parametrize these galactic and ex- 
tragalactic foregrounds with an eleven-parameter model. 
Five parameters describe the component amplitudes at 
i = 3000: three point source amplitudes in the 150 GHz 
band (clustering and Poisson for DSFGs and Poisson for 
radio) and a galactic cirrus amplitude in each of the two 
observation bands. One parameter describes the correla- 
tion between the tSZ and clustered DSFG terms. Three 
parameters describe the effective mean spectral indices of 
each point source component, and the final two param- 
eters describe the distribution about the mean spectral 
index for the Poisson radio source and DSFG popula- 
tions. 

The frequency dependencies of many of these fore- 
grounds are naturally discussed and modeled in units of 
flux density (Jy) rather than CMB temperature units. 
This is because the flux density of these foregrounds can 
be modeled as a power law in frequency, {Si, (x v°'). How- 
ever, in order to model these foregrounds in the power 
spectrum, we need to determine the ratio of powers in 
CMB temperature units. To calculate the ratio of power 
in the Vi cross Vj cross-spectrum to the power at (the 
arbitrary frequency) in units of CMB temperature 
squared, we multiply the ratio of flux densities by 



dB I dB I 
dT I t^o dT I '^0 
dB I dB I 
dT dT l^a 



(12) 



where B is the CMB blackbody specific intensity evalu- 
ated at TcmBj and and Vj are the effective frequencies 
of the SPT bands. Note that Vi may equal Vj . The ef- 
fective frequencies for each foreground calculated for the 
measured SPT bandpass (assu ming a nominal frequency 
dependence) are presented in ^6.4 



Poisson distributed point sources produce a constant 
Ci (thus, oc P). A population of dim (< 1 mJy) 
DSFGs is expected to contribute the bulk of the millime- 
ter wavelength point source flux. This expectation is con - 



sistent with recent observations (Lagache et al. 



Viero et al.|(|2009[), HIO, [Dunkley et al. (201017 



gOOT^, 
'nermal 

emission from dust grains heated to tens of Kelvin by star 
light leads to a large effective spectral index {Sv a J/"), 
wher e a > 2 between 15 and 220 GHz, for two reasons 
(e.g . [Dunne et alT 2000). First, these frequencies are in 
the Kayleigh- Jeans tail of the blackbody emission for the 
majority of the emitting DSFGs. Second, the emissivity 
of the dust grains is expected to increase with frequency 
for dust grains sm aller than the photon wavelength (2 
mm for 150 GHz) ( [Draine fc Lee[|1984[ [GOTdm[[T995| . 

Assuming that each galaxy has a spectral index a 
drawn from a normal distribution with mean a„ and vari- 



ance cTp, then following IMillea et al. ([2010|), the Poisson 
D5FU 



power spectrum of the Dsi' Gs can be written as 



„p Qp+0.51n(j;„j,„2)°"p 



3000 



(13) 



where r]i^i,,^2 — (^1^2/1^0) is the ratio of the frequencies of 
the spectrum to the base frequency, i'fgoo is the ampli- 
tude of the Poisson DSFG power spectrum a.t £ — 3000 
and frequency i^q. 

In principle, the frequency scaling of the Poisson power 
depends on both ap and Cp. HIO com bine the dus t spec- 
tral energy distributions (SEDs) from Silva et al. ( 1998 1 
with an assumed redshift distribution ot galaxies to in- 
fer (Tp 0.3 (although they conservatively adopt a prior 
range of ap G [0.2,0.7]). However, for a = 3.6, the ra- 
tio of Poisson power at 220 to 150 GHz only increases 
by 2% for ap = 0.3 versus ap = 0, which is small com- 
pared to the constraint provided by the current data. 
The complete model includes six parameters describing 
three measured Poisson powers in two bands. In order 
to reduce this to a reasonable number of parameters, we 
set ap — in the baseline model. The equivalent Poisson 
radio galaxy spectral variance, ar, has a negligible effect 
on the radio power at 220 GHz and is also set to zero. 

6.3.2. Clustered DSFGs 

Because DSFGs trace the mass distribution, they are 
spatially correlated. In addition to the Poisson term, this 
leads to a 'clustered' term in the power spectrum of these 
sources. 



(14) 



3000 



6.3.1. Poisson DSFGs 



DggQQ is the amplitude of the clustering contribution to 
the DSFG power spectrum a,t £ — 3000 at frequency i^o- 
The model for the clustering contribution to the DSFG 
power spectrum is which we have divided by 'I'looo 
to make a unitless template which is normalized at ^ = 
3000. 

We consider two different shapes for the angular power 
spectrum due to the clustering of the DSFGs, shown in 
Figure [5j The first is the shape of the fiducial model used 
in Hall et al. (2010). This spectrum was calculated as- 
suming that light is a biased tracer of mass fluctuations 
calculated in linear perturbation theory. Further, the 
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model assumes all galaxies have the same SED, that of a 
grey body, and the redshift distribution of the luminosity 
density was given by a parametrized form, with param- 
eters adjusted to fit SPT, BLAST, and Spitzer power 
spectrum measurements. Here we refer to this power 
spectrum shape as the "linear theory" template. 

On the physical scales of relevance, galaxy correla- 
tion functions are surprisingly well-approximated by a 
power law. The second shape we consider is thus a 
phenomenologically-motivated power law with index cho- 
sen to match the correlation p roperties observed for 
Lyman-break galaxie s at z ~ 3 (Giavalisco et al. 1998[ 
Scott fc White|1999| ). We adopt this "power law" model 
in the baseline model, and use a template of the form 
cx £'^'^ in all cases except as noted. 

We also consider a linear combination of the linear- 
theory shape and the power-law shape. Such a combi- 
nation can be thought of as a rough approximation to 
the two-halo and one-halo terms of a halo model. The 
present data set has limited power to constrain the ratio 
of the two terms; however, we do find a slight preference 
for the pure power-law model o ver th e pu re line ar-theory 
model, as discussed in sections |7.3.2| and |7.3.4| 

The spectrum of the clustered component should be 
closely related to the spectrum of the Poisson component 
since they arise from the same sources (HIO). However, 
the clustered DSFG term is expected to be more signifi- 
cant at higher redshifts, where the rest-frame frequency 
of the SPT observing bands is shifted to higher frequen- 
cies and out of the Raleigh- Jeans approximation for a 
larger fraction of the galaxies. One might therefore an- 
ticipate a slightly lower spectral index for the clustered 
term. 

We assume a single spectral index for the Poisson and 
clustered DSFG components in the baseline model {uc = 
Qfp), and explore the impact of allowing the two spectral 
indices to vary independently in |7.3.1| In general, we 
require ac = ap, except where otherwise noted. 

6.3.3. tSZ-DSFG correlation 

The tSZ signal from galaxy clusters and the clustered 
DSFG signal are both biased tracers of the same dark 
matter distribution, although the populations have dif- 
ferent redshift distributions. Thus, it is natural to ex- 
pect these two signals to be partially correlated and this 
correlation may evolve with redshift and halo mass. At 
frequencies below the tSZ null, this should result in spa- 
tially anti-correlated power. At 150 GHz, some models 
that attempt to associate emission with individual clus- 
ter member galaxies predict anti-correlations of up to 
tens of percent (SIO). LIO argue against this being sig- 
nificant in the SPT 150 GHz band, because cluster mem- 
bers observed at low redshift have significantly less dust 
emission than field galaxies. It is, of course, possible for 
this emission to be significant at higher redshifts where 
much of the tSZ power originates. Current observations 
place on ly w eak constraints on this correlation. (As we 
show in §7.4[ the data presented here are consistent with 
an anti-correlation coefficient ranging from zero to tens 
of percent.) To the extent that the Poisson and clustered 
power spectra have the same sp ectra l index, the differ- 
enced band powers discussed in ^7.5 (as well as in LIO) 
are insensitive to this correlation. 

A correlation, "f{t), between the tSZ and clustered 



DSFG signals introduces an additional term in the 
model, 



D 



tSZ-DSFG 



(15) 



We will assume no tS Z-DS FG correlation in all fits unless 
otherwise noted. In |7.4[ we will investigate constraints 
on the degree of correlation and the impact on other 
parameters from allowing a DSFG-tSZ correlation. 

6.3.4. Radio galaxies 

The brightest point sources in the SPT maps are coin- 
cident with known radio sources and h ave spectral indices 



consistent with synchrotron emission (Vieira et al. 2010 



hereafter VIO). Much like the Poisson Uyii G term, the 
model for the Poisson radio term can be written as 



Dl 



„r , ar+0.51ii(i7„^,„2)cr^ 



3000 



(16) 



The variables have the same meanings as for the treat- 
ment of the DSFG Poisson term. Clustering of the radio 
galaxies is negligible. 

In contrast to DSFGs, for which the number counts 
climb steeply towards lower flux, S^dN/dS for syn- 
chrotron sources is fairly constant. As a result, the 
residual radio source power (-D3000) ™^ the map depends 
approximately linearly on the flux above which discrete 
sources are masked. For instance, if all radio sources 
above a signal to noise threshold of 5g (6.4 m Jy) at 



150 GHz are masked, the de Zotti et al. ( |2005[ ) source 
count model predicts a residual radio power ot i'Jooo — 
1.28 /iK^. For a 10 cr threshold, the residual radio power 
would be i:»?;oop ^'2..b9^lK^. 

The de Zotti model is based on radio surveys such as 
NVSS extrapolated to 150 GHz using a set spectral in- 
dex for each subpopulation of radio sources. Extrapo- 
lating over such a large frequency range obviously intro- 
duces significant uncertainties; however, we do have an 
important cross-check on the modeling from the num- 
ber counts of bright radio sources detected in the SPT 
or ACT maps (VIO, [Marriage et~aI1 ( [20To| ). The de 
Zotti number counts appear to over-predict the number 
of very high flux radio sources reported by both experi- 
ments. We find that the de Zotti model should be mul- 
tiplied by 0.67 ± 0.14 to match the number of > 20mJy 
radio sources in the SPT source catalog (VIO). There 
are relatively few radio so urces at these fiuxes, and VIO 
and Marriage et al. (2010) have partially overlapping sky 
coverage. 'I'he counting statistics on these relatively rare 
objects will be improved significantly with the full SPT 
survey. The discrepancy disappears for lower flux radio 
sources and we find that, for 6.4 mJy < 6*150 < 20mJy 
150 GHz, the best fit normalization of the de Zotti model 
to fit the radio sources in the VIO catalog is 1.05 ± 0.14. 
Throughout this work, we use the de Zotti model predic- 
tion with a 15% uncertainty as a prior for residual power 
from sources below 6.4 mJy. Relaxing this ra dio prior 
by expandin g the 1 a amplitude prior about the |de Zotti] 
et al. (2005) model from 15% to 50% has no significant 



impact on resulting parameter constraints. 

We analyze the two-frequency VIO catalog to deter- 
mine the radio population mean spectral index (cxr) and 
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scatter in spectral indices (ar) at these frequencies. For 
each source, we take the full (non-Gaussian) likelihood 
distribution for the source spectral indices from VIO. We 
combine all sources above 5cr at 150 GHz with a prob- 
ability of having a synchrotron-like spectrum of at least 
50% in order to calculate the likelihoods of a given 
and ar- (We note that the results are unaffected by the 
specific choice of threshold probability for being a syn- 
chrotron source.) The best-fit mean spectral index is 
ar = —0.53 ± 0.09 with a 2 cr upper limit on the intrinsic 
scatter of ar < 0.3. With the present power spectrum 
data, both this uncertainty and intrinsic scatter are un- 
detectable. We have checked this assumption by allowing 
the spectral index to vary with a uniform prior between 
ar € [—0.66, —0.21] instead of fixing it to the best-fit 
value of -0.53 for several specific cases, and found no sig- 
nificant difference in the resulting parameter constraints. 
We therefore assume = —0.53 and = for all cases 
considered. 

6.3.5. Galactic cirrus 

In addition to the foregrounds discussed above, we in- 
clude a fixed galactic cirrus term parametrized as 



^£,1/1,1/2 ^ ^3000 



3000 



-1.5 



(17) 



HIO measured the Galactic cirrus contribution in the 
5^30™ fie ld by cross-correlat i ng the SPT map and model 
eight of [Finkbeiner et al.| ( |1999[ hereafter FDS) for 
the same held. We extend that analysis to the com- 
bined data set by scaling the 5*^30™ cirrus estimate 
by the ratio of the auto-correlation power in the FDS 
map of each field. We find the 23'^30'" field has 15% 
of the cirrus power measured in the 5*^30™ field. As 
the fields are combined with equal weight, we expect 
the effective cirrus amplitude in the combined band- 
powers to be 0.575 times the 5*^30™ amplitudes re- 
ported above. Based on this analysis, we fix the am- 
plitude of the cirrus term, cx i^^'^, in all param- 



eter fits to be {D 



cir,150,150 
3000 



p.cir,150,220 
-'^3000 



D 



cir,220,220 
3000 



) = 



(0.35, 0.86, 2.11) /iK^ respectively. Note that the cirrus 
power has a dust-like spectrum and is effecti vely r emoved 
in the differenced bandpowers discussed in j ]7.5| We as- 
sume there is zero residual cirrus power in the DSFG- 
subtracted bandpower analyses. Although we include 
this cirrus prior in the multi-frequency fitting, it has no 
significant impact on results; marginalizing over the am- 
plitudes in the two bands or setting them to either ex- 
treme (zero or the high value in the 5'^30™ field) does not 
change the resulting parameter values. 

6.4. Effective frequencies of the SPT bands 

Throughout this work, we refer to the SPT band center 
frequencies as 150 and 220 GHz. The data are calibrated 
to CMB temperature units, hence the effective frequency 
is irrelevant for a CMB-like source. However, the effec- 
tive band center will depend (weakly) on the source spec- 
trum for other sources. The actual SPT spectral band- 
passes are measured by Fourier transform spectroscopy, 
and the effective band center frequency is calculated for 
each source, assuming a nominal frequency dependence. 
For an q; = —0.5 (radio-like) source spectrum, we find 



band centers of 151.7 and 219.3 GHz. For an a = 3.8 
(dust-like) source spectrum, we find band centers of 154.2 
and 221.4 GHz. The band centers in both cases drop by 
^0.2 GHz if we reduce a by 0.3. We use these radio- 
like and dust-like band centers to calculate the frequency 
scaling between the SPT bands for the radio and DSFG 
terms. For a tSZ spectrum, we find band centers of 153.0 
and 220.2 GHz. The 220 band center is effectively at the 
tSZ null; the ratio of tSZ power at 220 to 150 GHz is 
0.2%. 

7. PARAMETER RESULTS 
7.1. Baseline model results 

As discussed in fj6] the baseline secondary anisotropy 
and foreground model contains four free parameters (be- 
yond the six ACDM parameters): the tSZ amplitude, 
the Poisson DSFG amplitude, the clustered DSFG am- 
plitude, and the DSFG spectral index. For the tSZ com- 
ponent, we use the SIO tSZ model scaled by a single 
amplitude parameter. For the clustered DSFG compo- 
nent, we use the power-law model presented in j ]6.3.2[ 
again scaled by an amplitude parameter. As discussed 
in j ]7.3.4[ the data show a slight preference for the power- 
law clustered point source shape model compared to the 
linear-theory model that was adopted in LIO. We require 
that both the clustered and Poisson DSFG components 
have the same spectral index, which we model as the free 
parameter ap, with no intrinsic scatter (tjp = 0). Table |3] 
shows the improvement in the quality of the fits with the 
sequential introduction of free parameters to the origi- 
nal ACDM primary CMB model. Beyond this baseline 
model, adding additional free parameters does not im- 
prove the quality of the fits. 

The additional parameters of the model discussed in 
f|6]arc all fixed or constrained by priors. The kSZ power 
is fixed to the D^qqq — 2.05 /iK^ expected for t he ho - 
mogeneous reionization model of Sehgal et al. (2010). 
Radio galaxies are included with a fixed spectral index 
ar — —0.53 and zero intrinsic scatter and an amplitude 
-^3000 = 1.28 /iK^ with a 15% uncertainty. Here, and 
throughout this work, we also i nclude a constant Galac- 
tic cirrus power as described in { 6.3.5 In addition to as- 
trophysical model components, for each band we include 
a term to adjust the overall temperature calibration and 
three terms to parametrize the uncertainty in the mea- 
surement of the beam function. This baseline model pro- 
vides the most concise interpretation of the data. As will 
be demonstrated in subsequent sections, relaxing priors 
on the fixed or tightly constrained parameters in this 
model has minimal impact on the measured tSZ power. 

Parameter constraints from this baseline model are 
listed in Table[2]and FigurejT] For the tSZ power, we find 
^30(X) ~ ^ ^-^ /iK ^, which is cons i stent with previous 
measurements (LIO, Dunkley et al. (2010l) and signifi- 
cantly lower than predicted by the SIO model. We detect 
both the Poisson and clustered DSFG power with high 
significance, with amplitudes of -Dfooo ~ '^■^ ^ ^.6 ^K^ 
and -D3000 — 6.1 ± 0.8 /iK^ respectively. The shared 
spectral index of these components is found to be a = 
3.58±0.09, in line wi th both theoretical exp ectations and 
previous work (HIO, Dunkley et al. (2010)). The ampli- 
tude of the Poisson term is significantly lower than what 
was measured in HIO; however, as we show in the follow- 
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TABLE 2 

Parameters from multi-frequency fits 



model 




nkSZ 




/ AAA 




etc 


3000 


r,C',220 

3000 


Baseline model 


3.5 ± 1.0 


[2.0] 


7.4 ±0.6 


6.1 ±0.8 


3.58 ±0.09 




70 ±4 


57 ±8 


Patchy kSZ 


2.9 ± 1.0 


[3.3] 


7.4 ±0.6 


5.7 ±0.8 


3.62 ±0.09 




72 ±4 


55 ±8 


Free kSZ 


3.2 ± 1.3 


2.4 ± 2.0 


7.4 ±0.6 


5.9 ± 1.0 


3.59 ±0.11 




71 ±5 


57 ±9 


Free ac 


3.4 ± 1.0 


[2.0] 


7.9 ± 1.0 


5.4 ± 1.4 


3.46 ±0.21 


3.79 ±0.37 


69 ±5 


59 ±8 


Linear-theory clustering 


3.5 ± 1.0 


[2.0] 


9.3 ±0.5 


4.9 ±0.7 


3.59 ±0.09 




89 ±5 


47 ±7 


Free ac &i linear theory 


3.4 ± 1.0 


[2.0] 


9.2 ±0.7 


5.3 ± 1.4 


3.62 ±0.14 


3.47 ± 0.34 


90 ±5 


46 ±7 


tSZ-DSFG correlation 


3.4 ± 1.0 


[2.0] 


9.3 ±0.7 


4.9 ±0.8 


3.60 ±0.13 




89 ±5 


46 ±7 


tSZ-DSFG cor, free kSZ 


1.8± 1.4 


5.5 ±3.0 


9.1 ±0.7 


4.4 ±0.9 


3.62 ±0.14 




89 ±5 


43 ±8 



Priors and Templates 

DG clustered 



model 


pars 


Free ac 


Free D^^o 


tSZ model 


kSZ model 


model 


-2Aln/: 


Baseline model 


10 


No 


No 


SIO 


SIO 


power law 




Patchy kSZ 


10 


No 


No 


SIO 


Patchy 


power law 


1 


Free kSZ 


11 


No 


Yes 


SIO 


SIO 


power law 





Free ac 


11 


Yes 


No 


SIO 


SIO 


power law 





Linear-theory clustering 


10 


No 


No 


SIO 


SIO 


linear theory 


4 


Free ac & linear theory 


11 


Yes 


No 


SIO 


SIO 


linear theory 


4 


tSZ-DSFG cor 


11 


No 


No 


SIO 


SIO 


linear theory 


4 


tSZ-DSFG cor, free kSZ 


12 


No 


Yes 


SIO 


SIO 


linear theory 


4 



Note. — Parameter constraints from multi-frequency fits to the models described in ^ The amplitudes are Vi at £ = 3000 
for the 150 GHz band, in units of /iK^. In all cases, the error specified is one half of tae 68% probability width for a given 
parameter after marginalizing over other parameters as described in the text. Where Qc is unspecified, both the clustered and 
Poisson DSFG components share an identical spectral index given by Up. The majority of the chains use a fixed kSZ spectrum; 
in these cases, the kSZ value appears in brackets. In addition to the model parameters described in ij6| we include two derived 

P 220 C 220 

parameters, -D30Q0 and -D3000 > which are the amplitude of the Poisson and clustered DSFG component in the SPT 220 GHz 
band, respectively. We also summarize the free parameters, templates, and relative goodness of fit for each case. 



ing sections, this discrepancy is eliminated by adopting 
the linear-theory clustered DSFG template used in that 
work. 

We explore the effect of altering this baseline model 
by replacing the tSZ template with the three alterna- 
tives discussed in §6.2.11 the Shaw, Trac, and Battaglia 
models. As anticipated, we find no significant difference 
between the models for any of the measured parameters. 
The resulting tSZ powers are listed in Table |6] they are 
statistically identical. The other model parameters also 
vary by less than 1% for all four tSZ model templates. 

7.2. kSZ variants 

The baseline kSZ model is the homogeneous reioniza- 
tion model presented in SIO. In this section, two alterna- 
tive kSZ cases are considered. First, for the same set of 
model parameters, we change the (fixed) kSZ power spec- 
trum from the baseline homogeneous reionization model 
to a template with additional power from patchy reion- 
ization. In the second case, we allow the kSZ amplitude 
to vary in order to explore the tSZ-kSZ degeneracy and 
to set upper limits on the SZ terms. 

We explore the effect of replacing the homogeneous 
kS Z temp late with the patchy kSZ template described 
The patchy kSZ model increases the power 
3.25 /iK^, and has a 



in ^SX^ 
at i 



3000 by 50% to L>f^of 



slightly modified angular scale dependence with addi- 
tional large-scale power. The likelihood of the best-fit 
model is nearly unchanged from the homogeneous kSZ 
case. The increase in kSZ power by 1.2 /xK^ leads to a 



decrease in tSZ power by 0.6 ^K^, as one would expect 
given the near-degeneracy along the tSZ + 0.5 kSZ axis 
which we discuss below. The kSZ spectrum shape differ- 
ence has no impact on the results. 

As an alternative to the fixed kSZ models discussed 
above, we also allow the kSZ amplitude to vary, adding 
one additional parameter to MCMC chains. With two 
frequencies, the current SPT data does not allow the 
separation of the tSZ, kSZ, and clustered point source 
components. The chains show a strong degeneracy be- 
tween the components, as shown in Figure IS] The best- 
constrained eigenvector in the kSZ/tSZ plane is approx- 
imately tSZ + 0.5 kSZ, very similar to the proportion- 
ality used in the differenced bandpowers in f 7.5 For 
this linear combination, we find D^qqq + 0.5^|qoq = 
4.5 ± 1.0 ^K^. We note that the data show some prefer- 
ence for a positive tSZ power, at just under 2a. We tested 
whether this preference is due to the detailed shapes of 
the two templates by swapping the tSZ and kSZ tem- 
plates. This model swap had no significant impact on 
the 2d likelihood surface, indicating that the preference 
for positive tSZ power is driven by the unique frequency 
signature of the tSZ. When marginalizing over the tSZ 
amplitude and other parameters, we find the upper limit 
on the kSZ power to be Df^oo < 6-5 fiK^ at 95% confi- 
dence. The corresponding upper limit on tSZ power is 
^30(fo < 5-3 /iK^ at 95% confidence. 



7.3. Clustered DSFG extensions to the baseline model 
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3) 



tsz 

3000 



TABLE 3 

Ax^ WITH THE ADDITION OF MODEL COMPONENTS 



Fig. 6. — The solid black bline shows the ID HkeUhood surface for £'3Q(X) from the multi-frequency analysis presented in [7.1 . The 
baseline model with the SIO homogeneous kSZ template has been assumed. Also plotted are the equivalent results from LIO (dot-dash 
line) and the DSFG-su btracted analysis in iffisl of this work (dashed line). In addition, the median value and 68% confidence region for 
the tSZ amplitude from |Dunkley et al.||20T(T| ) is shown by the diamond and thin horizontal line. 

The baseline model contains DSFG sources that gen- 
erate both Poisson and clustering contributions to the 
power spectrum. The spectral index of the Poisson power 
and clustering power is fixed to be the same. In this 
section, we examine two changes to this baseline model. 
First, we allow the clustered DSFG spectral index to vary 
independently, and second, we change the clustered tem- 
plate from a power-law to linear-theory model. We ex- 
plore all four permutations of these two Boolean cases 
(one of the four is the baseline case). One dimensional 
likelihood curves for each case are shown in Figure |9] 
and discussed below. We also consider a mixture of the 
power-law and linear-theory templates, each with its own 
variable amplitude, again fo r a to tal of five free parame- 
ters. As we will discuss in S 7.3.4 the SPT data slightly 
prefer the power-law template. 



Variable 


dof 


-2Aln£ 


point sources (Poisson) 


2 


-7821 


point sources (including clustering) 


1 


-104 


tsz 


1 


-25 


kSZ 





-1 



Note. — The addition of DSFG and SZ terms in the 
model significantly improves the quality of the maximum like- 
lihood fit. From top to bottom, the improvement in the fit 
as each term is added to a model including the primary 
CMB anisotropy and all the terms above it, e.g., the tSZ 
row shows the improvement from adding the tSZ effect with 
a free amplitude to a model including the CMB, and point 
source power, including the contribution from clustering of 
the DSFGs. The second column (dof), lists the number of de- 
grees of freedom in each component. The free Poisson point 
source parameters are the amplitude and spectral index of the 
DSFGs. Note that the third Poisson component (the radio 
power) is constrained by an external prior rather than the 
data and we have marginalized over the allowed radio power. 
The free clustered point source parameter is the amplitude of 
the power-law clustered template. The spectral index of the 
DSFG clustering term is fixed to that of the Poisson term. 
The last two rows mark the addition of the SIO tSZ template 
with a free amplitude, and the addition of the SIO homoge- 
neous kSZ template with a fixed amplitude. We do not find a 
significant improvement in the likelihood with any extension 
to this model, including allowing the kSZ amplitude to vary, 
or allowing independent spectral indices for the two DSFG 
components. 



As discussed in i 6.3.1 



7.3.1. Free 

the Poisson and clustered DSFG 



amplitudes are sensitive to different weightings of the 
DSFGs with redshift and brightness and thus we ex- 
pect these components t o have sligh t ly off set spectral 
indices. For instance, the Righi et al. (2008) model pre- 
dicts a^ — etc — 0.2. To study this oftset, we allow the 
two spectral indices to vary independently. The result- 
ing spectral indices for the clustered and Poisson DSFG 
terms separate by approximately one cr; howeve r, this 
shift i s in t he opposite direction from what the |Righi| 
et al. (2008) model predicts. With the addition of this 
additional tree parameter, the uncertainty on the Poisson 
and clustered amplitudes and spectral indices increase 
significantly. At 220 GHz, both DSFG powers are essen- 
tially unchanged, while the 150 GHz Poisson and clus- 
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Fig. 7. — ID likelihood plots for model parameters for three different kSZ treatments. The black, solid line is the baseline (fixed 
amplitude, SIO homogeneous model) kSZ model. The blue dotted line is the patchy kSZ model with fixed amplitude. The red, dashed 
line is the model where the amplitude of the (default) kSZ template is free to vary. All powers are measured at £ = 3000 and 150 GHz. 
The top left panel shows constraints on the tSZ power spectrum. The preferred tSZ power drops as the assumed kSZ power increases; this 
trade-off can be seen in more detail in Figure Is] The top right panel shows the likelihood curves for the Poisson DSFG power, and the 
bottom left panel shows the same for the clustered DSFG power. The bottom right panel shows the Poisson DSFG spectral index. There is 
a slight decrease in clustered DSFG power with the increased kSZ power expected from patchy reionization; this is accompanied by small 
shift to higher spectral indices, preserving the higher-frequency bandpowers. 



tered amplitudes slightly increase and decrease, respec- 
tively. 

7.3.2. Linear-theory model 

We adopt the HIO model as the linear-theory clus- 
tered DSFG model. The angular scale dependence of 
this model is very similar to the tSZ template and nearly 
flat in Dg. Swapping between the power-law and linear- 
theory models does not change the measured tSZ power. 
However, it does shift some power between the Poisson 
and clustered DSFG terms. The angular scale depen- 
dence of the power-law model is intermediate between 
the linear-theory model and Poisson terms. As a result, 
the flatter, linear-theory model leads to more Poisson and 
less clustered power. At £ = 3000 and 150 GHz, the Pois- 
son DSFG power increases by 3 cr to i^fooo — 9.3±0.5 fiK^ 
from 7.4 ± 0.6 ^K^. The clustered DSFG power de- 
creases correspondingly by 1 cr to -D3000 = 4-9 i 0-7 mK^- 
Note that the DSFG spectral index has dropped com- 
pared to what was measured in HIO from 3.86 ± 0.23 
to Up = 3.59 ± 0.09. The change in spectral index is 
driven by the 220 GHz power: the Poisson power is 2 cr 
lower and the clustered power is 1 cr lower than what was 
measured in HIO. 

7.3.3. Free olc & linear-theory model 



Finally, we consider the case where we allow ttc as a free 
parameter with the linear-theory clustered DSFG tem- 
plate. The resulting chains include features from both 
the free-ctc and linear-theory cases: the Poisson power 
increases and tSZ power decreases, while uncertainties 
on the clustered amplitude and all spectral indices in- 
crease. The spectral indices of the clustered and Poisson 
DSFG terms do not separate as much as they do in the 
power-law, free-ctc chain. 

7.3.4. Discriminating between the power-law and 
linear-theory models 

In principle, the likelihood of a fit to the data can be 
used to distinguish between the power-law and linear- 
theory clustered DSFG templates described above. We 
explore this in two ways. First we consider the improve- 
ment in the best-fit log-likelihood between the models 
which use each template. The power-law model produces 
a marginally better fit to the data, with Aln£ — —2 when 
compared to the linear-theory model. 

In addition to the choice of one template or the other, 
we extend the baseline model to include a second DSFG 
clustering component. The two DSFG clustered ampli- 
tudes, corresponding to the power-law and the linear- 
theory templates, are allowed to vary independently. The 
spectral index for all three (Poisson and two clustered) 
DSFG components are tied. The two-dimensional like- 
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Fig. 8. — 2D likelihood contours for r'Hoo 



vs. fjjoQO' 68%, 95% and 99.7% contours are marked with blue, orange, and gray 

respectively. Both tSZ and kSZ model templates are from SIO. As in the differenced spectrum analysis, the data primarily constrain 
a linear combination Dt^Q^^ + O.SDj^qq = 4.5 it 1.0 fiK^. The dashed line shows this axis of near degeneracy. Due to this degeneracy, 
neither parameter is significantly detected above zero. However, 95% CL upper limits of 5.3 and 6.5 /iK^ can be placed on the tSZ and kSZ 
respectively. 



lihood surface for the amplitudes of each template is 
shown in Figure 10 The data show a slight preference for 
the power-law template over the linear-theory template; 
however, we can not rule out the linear-theory model or 
a linear combination of the two with any confidence. We 
can expect this constraint on DSFG clustering shape to 
improve with future data. 



7.4. tSZ-DSFG 



slation 



As discussed in i 6.3.3 there are reasons to expect that 
the signal from the tSZ and the clustered DSFGs may be 
correlated. In order to study that correlation, we extend 
the model to include the correlation coefficient of these 



components, detailed in ^6.3.3 We assume the cor- 
relation coefficient is constant in i, thus reducing j{£) 
to 7. We would expect an anticorrclation at frequencies 
below the SZ null; however, we do not require that the 
correlation coefficient is negative. Instead, we set a flat 
prior on the correlation 7, and allow it to range from -1 
to 1. 

We consider the impact of introducing correlations be- 
tween the tSZ and DSFGs on two chains - both allow 
the correlation factor to vary and one also allows the kSZ 
amplitude to vary. In both chains, we assume the linear- 
theory clustered DSFG template, because its shape is 
similar to the SIO tSZ template. 

The results for this model without the possibility of a 
tSZ-DSFG correlation have been presented in [ 7.3.2 The 
marginalized Id likelihood function for the tSZ-DSFG 
correlation factor is shown in Figure |11[ In the fixed 
kSZ case, the data are consistent with zero correlation: 
7 = 0.0 ± 0.2. When the kSZ power is allowed to vary, 
the data remain consistent with zero correlation, how- 
ever, large anti-correlations are also possible. The 95% 



CL upper limit on D^qqq allowing for a tSZ-DSFG cor- 
relation is unchanged from the equivalent chain with the 
correlation coefiicient set to zero; marginalizing over a 
non-zero tSZ-DSFG correlation factor does not allow ad- 
ditional tSZ power. As Figure [T2| shows, this is because 
kSZ power increases for larger negative correlation in 
order to preserve the observed power in the 150 x 220 
bandpowers. Naively, one might expect the best fit tSZ 
power to increase when the possibility of a tSZ canceling 
correlation is taken into account; however, in practice, a 
negative correlation leads to less tSZ power in the best fit 
model. The anti-correlation reduces the predicted power 
in the 150 x 220 spectrum, and kSZ power is increased to 
compensate. The sum of kSZ and tSZ is well constrained 
by the data, so more kSZ power leads to less tSZ power. 
When marginalizing over the tSZ-DSFG correlation, the 
95% CL on tSZ decreases slightly to DI^§q < 4.5 ^K^, 
while the kSZ upper limit increases to -Dfifoo < 9.7/iK^. 
From this, as well as the agreement between the multi- 
frequency fits and the DSFG-subtracted result presented 
in §7 5| we conclude that the addition of a tSZ-DSFG 
correlation cannot explain the discrepancy between the 
observed tSZ power and models that predict a signifi- 
cantly larger tSZ contribution. 

These results depend on the Poisson and clustered 
DSFG terms having similar spectral indices. We expect 
the upper limits to be substantially weakened if the two 
spectral indices are allowed to separate, but note that 
the required spectral index separation would need to be 
much larger than the ap — Uc ^0.2 difference expected 
from models (HIO). 

7.5. Comparing results from differenced spectra to 
multi-frequency fits 
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Fig. 9. — ID likelihood plots of model parameters for different DSFG clustering treatments. Clockwise from top left: the ID likelihood 
curves for the amplitudes of the tSZ effect, Poisson DSFG amplitude, Poisson DSFG spectral index and clustered DSFG amplitude. Four 
cases are plotted for each parameter: The black solid line indicates the baseline model, where ac = ccp and the power-law clustered 
template is used. The red dashed line indicates the case where Oc is a free parameter and the power-law clustered template is used. 
The blue dotted line indicates the case where Qc = cip and the linear-theory clustered template is used. The blue-green dot-dash 
line indicates the case where «c and ap are independent and the linear-theory clustered template is used. Using the power-law clustered 
template tends to reduce the power in the Poisson component compared to the linear-theory template. This is expected as the power-law 
shape has more power at high-£ than the linear-theory template and the Poisson component also increases with £ (Di oc i'^). Including an 
additional free parameter (etc) tends to increase the uncertainty on DSFG-related parameters but not on tSZ. 



The main result of this work is the set of multi- 
frequency bandpowers. However, we also follow the 
treatment in LIO and consider a linear combination of 
the two frequency maps designed to eliminate the fore- 
ground contribution from DSFGs. Specifically, this is the 
power spectrum of the map m constructed from the 150 
and 220 GHz maps {mi^Q and m22o respectively) accord- 
ing to 
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1 - X 



("1-150 -xm22o)- 



(18) 



In LIO, the subtraction coefficient x was selected by 
minimizing the Poisson point source power in the resul- 
tant bandpowers. We repeat this analysis for the new 
bandpowers and find x — 0.312 ± 0.06. This is consis- 
tent with the value oi x — 0.325 ± 0.08 reported in LIO. 
For complete details on the analysis, we refer the reader 
to that work. Given the consistency of these results, 
we have chosen to use x ~ 0.325 to create the 'DSFG- 
subtracted' bandpowers presented in Table |4] to enable 
the direct comparison of the two data sets. These two 
subtraction coefficients correspond to a spectral index 
for the DSFG population of ap = 3.7 and 3.6 respec- 
tively, consistent with the ap = 3.57 ± 0.09 determined 
in multi-frequency fits. We compare the LIO bandpowers 
with those presented in this work in Figure 13 and find 



excellent agreement between them. Note that the two 
data sets are not independent since LIO analyzed one of 
the two fields used in this work. 

This 'DSFG-subtracted' power spectrum is interesting 
for two reasons. First, it allows a direct comparison with 
the results of LIO. Second, in the DSFG-subtracted 
bandpowers, the residual correlated DSFG power will 
only be significant if the frequency spectrum of DSFGs 
in galaxy clusters is substantially different than the spec- 
trum of DSFGs in the field. Therefore, parameter fits 
to DSFG-subtracted bandpowers are independent of the 
assumed DSFG power spectrum shape or tSZ-DSFG cor- 
relation, and are a potentially interesting cross-check on 
the multi-frequency fits. However, in practice, it is worth 
noting that these fits constrain a similar linear combina- 
tion of tSZ and kSZ, roughly given by tSZ -I- 0.5 kSZ. 
Both methods therefore result in a similarly robust de- 
tection of tSZ power independent of tSZ-DSFG correla- 
tion. 

For this DSFG-subtracted case, we simplify the base- 
line model to include only one free parameter beyond 
ACDM: the tSZ power spectrum amplitude. We also in- 
clude a Poisson point source term with a strict prior, and 
a fixed kSZ power spectrum. The SIO model is used for 
both kSZ and tSZ. The Poisson prior is identical to the 
point source prior described in LIO with one exception: 
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TABLE 4 

DSFG-SUBTRACTED BANDPOWERS 



i range 


^cff 




G (/iK^) 


2001 - 


2200 


2056 


245.0 


9.6 


2201 - 


2400 


2273 


140.3 


6.9 


2401 - 


2600 


2471 


110.5 


6.0 


2601 - 


2800 


2674 


59.8 


5.1 


2801 - 


3000 


2892 


59.0 


5.4 


3001 - 


3400 


3185 


38.9 


3.9 


3401 - 


3800 


3581 


23.8 


4.3 


3801 - 


4200 


3994 


20.7 


5.4 


4201 - 


4600 


4402 


20.9 


6.7 


4601 - 


5000 


4789 


12.3 


7.9 


5001 - 


5900 


5449 


17.4 


7.3 


5901 - 


6800 


6360 


13.0 


11.8 


6801 - 


7700 


7257 


60.7 


18.4 


7701 - 


8600 


8161 


3.6 


28.6 


8601 - 


9500 


9063 


56.2 


44.8 



Note. — Band multipole range and weighted value £cff, bandpower D, and uncertainty a for the DSFG-subtracted auto- 
spectrum of the SPT fields. This is the power spectrum of a linear combination of the 150 and 220 GHz maps ((miso — 
0.325 m22o)/(l — 0.325)) designed to suppress the DSFG contribution. The quoted uncertainties include instrumental noise and 
the Gaussian sample variance of the primary CMB and the point source foregrounds. The sample variance of the SZ effect, 
beam uncertainty, and cali brat ion uncertainty are not included. Beam uncertainties are shown in Figure [T] and calibration 
uncertainties are quoted in §2.2| Point sources above 6.4 mjy at 150 GHz have been masked out in this analysis, eliminating the 
majority of radio source power. 
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Fig. 10.— 68% (blue), 95% (orange) and 99.7% (grey) likeli- 
hood contours for the amplitudes of clustered DSFG components 
for a model containing both power-law and linear-theory com- 
ponents with independently variable amplitudes. The power-law 
model is favored at the 1 u level; the data do not clearly discrimi- 
nate between the power-law shape, linear-theory shape, or a linear 
combination of the two. The choice of clustered DSFG model has 
no impact on the tSZ power and we adopt the power-law template 
as the baseline model. 

the radio source aniplitude prior has been tightened to 
match the 15% prior used in this work instead of the 
50% log-normal prior used in LIO. However, tests with 
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Fig. 11.— The ID likelihood curve for the tSZ-DSFG correlation 
factor, 7. The linear-theory clustered DSFG template and the SIO 
tSZ template are assumed. The solid line indicates the case where 
the kSZ amplitude is held fixed. The data favor no correlation al- 
beit with significant uncertainty. The dashed line denotes the 
likelihood curve when the kSZ amplitude is allowed to vary. The 
data remain consistent with zero correlation, but large negative 
correlations become possible. Essentially, the kSZ power increases 
to balance the lower predicted power due to anti-correlation. Be- 
cause the tSZ and kSZ amplitudes are highly correlated, additional 
kSZ power reduces the apparent tSZ power. 



relaxed priors show the results are insensitive to this ra- 
dio uncertainty. Except for this small change in the radio 
prior, and the switch to WMAP7 data, these differenced 
spectrum fits are identical to the parameter fitting used 
in LIO. 

As shown in Figure [6j the resulting DSFG-subtracted 
tSZ amplitude agrees with both the LIO value and the 
multi-frequency result presented here. We measure the 
tSZ power at 150 GHz to be Df^^^ = 3.5 ± 1.2 /.tK^, 
with a fixed (SIO homogeneous) kSZ component. For the 
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Fig. 12.— 2D likelihood contours for tSZ-DSFG correlation chains in which kSZ is allowed to vary. The 68%, 95% and 99.7% likelihood 
contours are marked with blue, orange and gray respectively. From left to right: tSZ power vs. correlation coefficient 7; tSZ power vs. kSZ 
power, marginalized over 7 (and foreground terms); and kSZ power vs. 7. 
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Fig. 13. — The SPT DSFG-subtracted bandpowers from this work (black circles) versus the DSFG-subtracted bandpowers presented 
in LIO (gray squares), an earlier reduction of one half of the 2008 data. The damping tail of the primary CMB anisotropy is apparent 
below I = 3000. Above I = 3000, there is a clear excess due to secondary anisotropics and residual point sources. The two datasets show 
excellent consistency. 



same template choices, LIO report I^Ioto = 3.2±1.6 /^K^, 
and the multi-frequency fit to the baseline model finds 
^30(X) — 3.5 ± 1.0 ^K^. This close agreement demon- 
strates the consistency of the two methods of parameter 
fitting, and confirms the results presented in LIO. 

8. CONSEQUENCES OF THE MEASURED tSZ SIGNAL 

We now focus on comparing the improved measure- 
ments of the tSZ power spectrum with the predictions of 
recent models and simulations. We then combine mea- 
surements of the primary CMB power spectrum with the 
tSZ signal to obtain joint constraints on a^. 

8.1. Thermal SZ amplitude 



LIO found the measured amplitude of the tSZ power 
spectrum to be a factor of 0.42 ± 0.21 lower than that 
predicted by the simulations of SIO (at their fiducial cos- 
mology). In this analysis, having doubled the survey 
area to 210 deg^, we find that the amplitude of the tSZ 
power spectrum at ^ = 3000 is L'gf^o = 3-5 ± 1.0 /xK^ at 
150 GHz. This is 0.47 ± 0.13 times the SIO template. 

The discrepancy suggests that either the SIO model 
overestimates the tSZ power, or the assumed set of 
cosmological parameters may be wrong. To calculate 
the tSZ power spectrum, SIO combined a semi-analytic 
model for intra-group and cluster gas with a cosmological 
N-body lightcone simulation charting the matter distri- 
bution out to high redshift. The gas model assumed by 
SIO was calibrated by comparing it with X-ray observa- 
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Fig. 14. — Poste rior distribution of the tSZ scaling parameter 
^tSZ (see Equation |19| and Table|6| for each of the tSZ templates. 
The black line shows the SIO preaiction. The red solid, dashed 
and dot-dashed lines show the results for the Shaw, Trac, and 
Battaglia models, respectively. As expected, the preferred value of 
Atsz increases for models which predict less tSZ power. Note that 
a value of Atsz = 1 would demonstrate consistency between the 
predicted and observed tSZ power. 

tions of low-redshift (z < 0.2) clusters. However, the ma- 
jority of tSZ power is contributed by groups or clusters 



at hi gher redshift (Komatsu & Seljak 2002 Trac et al. 



2010| . Hence, the yiU model may substantially over 
predict the contribution from higher redshift structures. 
Note that in §7.4[ we demonstrated that the lower-than- 
predicted tSZ power is unlikely to be due to a spatial 
correlation between DSFGs and the groups and clusters 
that dominate the t SZ sig nal. 

As described in §6.2. 1[ in addition to the SIO tem- 
plate we inve stigate resul t s obta i ned using th e temp lates 



I 

tE 



suggested by|Trac et al.| (|2010|), |Shaw et al.| ( |2010[ ) and 

tE 



2010 ((the I'rac, Shaw, and Hattaglia 



Battaglia et aT 

models). I'hese three alternative models all predict sub 
stantially less power aX £ = 3000 than the SIO template 
(assuming the same cosmological parameters). For each 
template we run a new chain substituting the SIO tem- 
plate with one of the three others. The rest of the base- 
line model remains as described in 36] 

Alternatively, the amplitude of the tSZ power spec- 
trum is extremely sensitive to erg, scaling approximately 
as D*/^ cx al~^{flbh)^. The exact dependence on erg de- 
pends on the relative contribution of groups and clusters 
as a function of mass and redshift to the power spectrum, 
and thus the a strophysic s of the intra-clus t er medium 
(iKomatsu &: Seljak 2002) |Shaw et al.]|2010[ [Trac et al.| 
2010p . All tour of the templates that we consider in tEis 
work were determined assuming erg = 0.8. The deficit 
of tSZ power measured in the SPT maps may therefore 
indicate a lower value of CTg . 

To correctly compare the measured tSZ power spec- 
trum amplitude with theoretical predictions, we must 
rescale the template power spectra so that they 
are consistent with the cosmological parameters at 
each point in the chains. Recall that the cosmo- 
logical parameters are derived solely from the joint 
WMAP7-|-QUaD-HACBAR-|-SPT constraints on the pri- 
mary CMB power spectrum and so are not directly in- 



fluenced by the amplitude of the measured tSZ power. 

We define a dimensionless tSZ scaling parameter, Atsz, 
where 

ntsz 

^tsz(..) = -^f^ . (19) 



(f>tSZ 



ntsz 

^3000 
^3000 1 



is 



the measured tSZ power at £ — 3000 and 
Ki) is the template power spectrum normalized to 
be consistent with the set of cosmological parameters 'k^' 
at step i in the chain. We obtain the preferred value of 
^tsz by combining the measurement of D^qqq with the 
cosmological parameter constraints provided by the pri- 
mary CMB measurements. 

To determine the cosmological scaling of each of the 
templates, we ad opt the analytic model described in 
Shaw e t aL|(2010 ). This model uses the halo mass func- 
tion ot I l inker et al. 



( 2008 1 to calculate the abundance 
of dark matter halos as a function of mass and redshift, 
combined with an analytic prescription (similar to that of 
Bode et al. | 2009[) to determine the SZ signal from groups 



and clusters. Snaw et al. (2010) demonstrated that when 



the same assumptions are made about the stellar mass 
content and feedback processes in groups and clusters, 
their model almost exactly reproduces the SIO template. 
For the Trac and Battaglia models, we modify the astro- 
physical parameters of the Shaw model so that it agrees 
with the other models at their fiducial cosmology. Vary- 
ing the input cosmological parameters thus allows us to 
estimate the cosmological scaling of the tSZ power spec- 
trum for any given model. 

As an indication of the relative cosmological scaling 
of each model a.t £ = 3000, in Table [5] we give the de- 
pendence of $tsz on erg, fib, and h around their fiducial 
values. We assume that the cosmological scaling takes 
the form. 



^tsz 

^3000 



oc 



f-V 

\0.8J 



0.044 



h 
071 



(20) 



Note that the values of a given in Table [5] for the SIO and 
Trac m odels agree well with those given in |Trac et "ah 
(20101. We find the templates depend only weakly on 
other cosmological parameters at this angular scale. 



TABLE 5 
tSZ Cosmological Scaling 



model 


a 




7 


Sehgal 


8.09 


2.50 


1.89 


Shaw 


8.34 


2.81 


1.73 


Trac 


8.34 


2.90 


1.81 


Battaglia 


8.33 


2.85 


1.71 



Note. — a, /3, and 7 describe the cosmological scaling 
of the tSZ power spectrum for each template, as defined 
in Equation [20] 



Figure [14] shows the posterior probability of Atsz ob- 
tained forThe SIO, Shaw, Trac, and Battaglia templates. 
The results are also given in Table [6] Recall that a value 
of Atsz = 1 would demonstrate that the model is con- 
sistent with the measured tSZ power, while Atsz less 
(greater) than 1 indicates that the model over- (under-) 
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estimates the tSZ power spectrum amplitude. All four 
models display a non-Gaussian distribution with tails ex- 
tending to higher values of ^tsz- Therefore, in Table [6] 
we give both symmetric 68% confidence limits for both 
Atsz and ln(y4tsz) as well as the asymmetric 95% confi- 
dence limits (in brackets) for A^sz- 

For the SIO template, the preferred value is Atsz = 
0.33 ± 0.15. The measured level of tSZ power is thus sig- 
nificantly less than that predicted by the SIO template. 
When the kSZ amplitude is also allowed to vary, we find 
a slight decrease in Atsz to 0.30 ± 0.16. Note that the 
constraint on ag obtained from the primary CMB power 
spectrum is 0.821 ±0.025, which is greater than the fidu- 
cial value of 0.8 assumed by all the templates}^ The 
values of Atsz are therefore less than what we would 
have obtained had we fixed the templates at their fidu- 
cial amplitudes. 

We find that the discrepancy remains, but is reduced, 
for the Shaw, Trac, and Battaglia models. The Trac 
and Shaw templates provide the best match to the mea- 
sured amplitude with Atsz = 0.59 ± 0.25. It is interest- 
ing to note that the total amount of tSZ power measured 
(-^30(X) Table ^ is essentially independent of the tem- 
plate used. This aemonstrates that the differences in the 
shapes of the templates do not influence the results. 

The templates that are most consistent with the ob- 
served SPT power at £ = 3000 (when combined with 
cosmological constraints derived from the primary CMB 
power spectrum) are those that predict the lowest tSZ 
signal. The SIO, Trac, and Shaw models adopt a similar 
approach for predicting the tSZ power spectrum. Each 
assumes that intra-cluster gas resides in hydrostatic equi- 
librium in the host dark matter halo, and include sim- 
ple prescriptions to account for star formation and AGN 
feedback. The Trac and Shaw models - which both pre- 
dict significantly less tSZ power - differ from the SIO 
model in two princip al ways. First, follow ing the observa- 
tional constraints of 
greater stellar mass 



Giodini et al. ( 2009 ) , they assume a 
raction m groups and clusters. This 
reduces the gas density and thus the magnitude of the SZ 
signal from these structures. Second, the Trac and Shaw 
models include a significant non-thermal contribution to 
the total gas pressure. Hydrodynamical simulations con- 
sistently observe significant levels of nonthermal pressure 
in the intra- cluster medium due to bulk gas motions and 



Rasia et al. 

2om 



2004 Kay et al. 



In terms of the 



turbulence (Ev rard' 1990 
2004||Dolag"eral. 2005 ; Lau et al 
tyz power spectrum, including a non-thermal pressure 
component reduces the thermal pressure required for the 
ICM to be in hydrostatic equilibrium, therefore reducing 
the predicted tSZ signal. 

The Battaglia template was produced using maps con- 
structed from a hydrodynamical simulation that included 
star formation and AGN feedback. While direct compar- 
isons with the mode l parameters of the oth er three tem- 



plates are difficult, Battaglia et al. (2010) found that. 



in simulations in which AGIN feedback was not included 
the tSZ power spectrum increased by 50% at the angular 
scales probed by SPT. If we assume the tSZ template 

^1 We note that, following the upgrade to RECFAST vl.5 and in 
the WMAP likelihood code to v4.1, the WMAP7-al one preferred 
value of erg is O.Sll^Q Qg^ as given in Appendix B of 
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Fig. 15. — 2D likelihood contours at 68% and 95% con fidence 
for us vs. the tSZ scaling parameter A^sz (see Eqn. |19| . Red, 
dashed contours show results having resampled the chains ap- 
plying a lognormal prior to Atsz centered at In(Atsz) = and 
of width o"in(Atsz) = 0.085. This prior is chosen to reflect the ex- 
pected sample variance for the observed sky area (see text). Black 
contours show the results without this prior. All results are for 
the baseline model (SIO tSZ model and fixed homogeneous kSZ). 
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Fig. 16. — erg constraints for the SIO (black line), Shaw (red, 
solid line), Trac (red, dashed line) and Battaglia (red, dot- 
dashed line) models. No model uncertainty has been included. 

produced from their simulations without AGN feedback 
we find Atsz — 0.27±0.10 and thus a significant increase 
in the disc repency between the o bserved and predicted 
tSZ power. Battaglia et al. (2010) demonstrate that the 
inclusion of AGN feedback reduces the thermal gas pres- 
sure in the central regions of clusters by inhibiting the 
inflow of gas. This results in flatter radial gas pressure 
profiles and a suppression of tSZ power on angular scales 

e > 2000. 

8.2. (Tg Constraints 

To obtain joint constraints on erg from both the pri- 
mary CMB anisotropics and the tSZ amplitude we must 
determine the probability of observing the measured 
value of Atsz for a given set of cosmological parameters. 
In the absence of cosmic variance - and assuming the 
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TABLE 6 
tSZ Constraints 



model 






In(Asz) 


0-8 


Sohgal 
Sehgal, free kSZ 
Shaw 
Trac 
Battaglia 


3.5 ± 1.0 
3.3 ± 1.3 

3.6 ± 1.0 
3.6 ± 1.0 
3.5 ±0.9 


0.33 ± 0.15 (Igig) 
0.30 ± 0.16 (IJJ jg) 
0.58 ± 0.27 

0.59 ± 0.25 (t;;.!?) 

0.45 ±0.20(^0.3?) 


-1.14 ± 0.45 
-1.27± 0.47 
-0.54 ± 0.40 
-0.53 ±0.45 
-0.79 ±0.44 


0.771 ±0.013 
0.772 ± 0.013 
0.799 ±0.014 
0.802 ±0.013 
0.786 ±0.014 



Note. — Measured tSZ power and inferred as const rain ts for different tSZ and kSZ model assumptions. Note that the 
posterior distribution of Atsz is non-Gaussian (see Figure 141. To account for this we give both the 68% and asymmetric 95% 
(in brackets) confidence intervals. We also give the median and 68% confidence interval for In(yltsz). For Dlow ^nd erg we quote 
only the 68% confidence interval. The erg constraints do not account for model uncertainty. All results assume the homogeneous 
reionization kSZ template given in SIO. In the 'free kSZ' case the amplitude of this template is aUowed to vary, otherwise it is 
held fixed at 2.05 /iK^. 



predicted tSZ template to be perfectly accurate - the 
constraint on erg would be given by its distribution at 
^tsz = 1 for each of the models. However, in practice, 
the expected sample variance of the tSZ signal is non- 
negligible. Based on the hydrodynamical simulations of 
Shaw et al. (2009), LIO assumed the sample variance at 
I = 30UU to be a lognormal distribution of mean and 
width crin(Atsz) ~ 0.12. The survey area analyzed here 
is twice that of LIO and thus we reduce the width of 
this distribution by a factor of v^- For each of the four 
tSZ templates, we construct new chains to include this 
prior by importance sampHng. We address the issue of 
theoretical uncertainty in the tSZ templates later in this 
section. 



The black contours in Figure 15 show the 68% and 
95% confidence limits in the cts — Atsz plane for the SIO 
model before applying the prior on Atsz- In this case, 
the constraint on erg is derived entirely from the primary 
CMB power and is not infiuenced by the measured tSZ 
signal. The combined WMAP7±ACBAR±QUaD+SPT 
constraints on the primary CMB yield erg = 0.821±0.025. 
The shapes of the contours reflect the dependence of Atsz 
on (jg; higher values of erg predict more tSZ power and 
require a lower tSZ scaling parameter, whereas lower val- 
ues of tJg suggest a tSZ amplitude closer to that predicted 
by the SIO template. 



The dashed red contours in Figure 15 represent the 
joint 68% and 95% confidence limits on erg and Atsz for 
the SIO template having resampled the chain to include 
the prior on Atsz- The ID e rg c onstraints for all four tem- 
plates are shown in Figure [16] and Table |6] For the SIO 
template the preferred value~ol erg is 0.771 ±0.013j^ The 
Trac, Shaw, and Battaglia models all prefer higher values 
of erg . The Trac template - the most consistent with the 
measured tSZ amplitude - gives erg = 0.802±0.013. This 
represents nearly a factor of two improvement in the sta- 
tistical constraint on erg compared to that obtained from 
the primary CMB alone. 

It is clear from Table [6] that the four templates pro- 
duce a spread in the preferred value of erg that is larger 

During the analysis for this work we discovered that the joint 
primary CMB-tSZ constraint on (jg in |Lueker et aL] | |2010[ l erro- 
neously assumed the SIO model was created using erg = u!77 in- 
stead of the correct value of trg = 0.8. The constraints on erg should 
therefore have be en 0.767± 0.018 (as suming no theory uncertainty, 
see Section 6.4 in |Lueker et al.|2010| and 0.778 ± 0.022 (assuming 
a 50% theory uncertamty ) . 



than the statistical uncertainty on any of the individual 
results. This spread reflects the range in amplitude of 
the tSZ templates (at a fixed value of erg) due to the dif- 
ferences in the thermal pressure profiles of groups and 
clusters in the underlying models and simulations. The 
dominant uncertainty on the tSZ-derived erg constraint is 
thus the astrophysical uncertainty on the predicted tSZ 
amplitude rather than the statistical error on the mea- 
surement of -Daooo- Accounting for this uncertainty will 
clearly degrade the constraint on erg. 

To produce a constraint on erg that accounts for the 
range in amplitude of the tSZ templates we must add 
an additional theory uncertainty to the prior on Atsz- 
Assuming the Trac template - which is the most consis- 
tent with the measured tSZ signal - we find that a 50% 
theory uncertainty on In(Atsz) is sufficient to encom- 
pass the other three templates. As the sample variance 
is insignificant compared to the theory uncertainty, the 
new prior on Atsz is effectively a lognormal distribution 



of width (Tin 



(Atsz) 



0.5. With this additional uncer- 



tainty, the revised constraint on erg for the Trac template 
is 0.812 ± 0.022. The astrophysical uncertainty clearly 
degrades the constraint significantly - it now represents 
only a small reduction of the statistical uncertainty of 
the primary CMB-only result. 

It is clear that, while the amplitude of the tSZ power 
spectrum is extremely sensitive to the value of erg, the 
constraints that can obtained on this parameter are lim- 
ited by the accuracy with which we can predict the 
tSZ power spectrum amplitude for a given cosmological 
model. The principal difficulty in producing precise pre- 
dictions is that groups and clusters spanning a wide range 
of mass and redshift contribute to the power spectrum 
at the angu l ar sca les being probed by SPT. For instance. 



Trac et al.| (2010) find that 50% of the power in their 
nonthermal2U template is due to objects at z > 0.65. 
Similarly, they flnd that 50% of the power comes from 
objects of mass less than 2 x 10^*h^^ Mq. To date, 
both these regimes have been poorly studied by targeted 
X-ray and SZ observations. It is therefore difficult to 
test current models and simulations of the tSZ power 
spectrum sufficiently well to produce precise predictions. 
X-ray and SZ observations of groups and high redshift 
clusters will provide constraints on models and simula- 
tions that should improve estimates of the tSZ power 
spectrum and thus allow tighter constraints on erg from 
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2010). Alterna- 



measurements of this signal ( Sun et al. 
tively, better constraints on cosmological parameters and 
improved measurements of the tSZ power spectrum will 
continue to shed new light on the relevant astrophysics 
in low mass and high redshift clusters. 

9. CONCLUSIONS 

We have presented the CMB temperature anisotropy 
power spectrum from the complete 210 deg^ of sky ob- 
served by SPT in the 2008 season. Bandpowers from 
observations at 150 and 220 GHz are shown in Table [TJ 
The bandpowers correspond to angular scales, I > 2000, 
where secondary anisotropics and foreground contribu- 
tions are significant. We also present bandpowers that 
have been spectrally differenced to remove the power 
contributed by dusty sources, following the method de- 
scribed in the first SPT power spectrum release (LIO). 

We perform multi-frequency model fits to the £ > 2000 
bandpowers to constrain secondary CMB anisotropics 
and foreground signals. The minimal extension to the 
six-parameter lensed ACDM primary CMB spectrum 
that we consider consists of four free parameters: the 
amplitudes of the tSZ spectrum, Poisson and clustered 
DSFG terms, and the spectral index of the DSFG 
sources. The addition of these four parameters improves 
the likelihood of the best-fit model by Aln£ = —3980. 
This baseline model also includes a fixed kSZ contribu- 
tion and a tight prior on radio sources. We explore the ef- 
fect of modifying this model by allowing additional model 
parameters to vary. In particular, we free the amplitude 
of the kSZ and the clustered DSFG spectral index. How- 
ever, we do not find that these additional free parameters 
improve the quality of fits to the data. 

With these data, the amplitudes of the two SZ com- 
ponents are largely degenerate. The data primarily con- 
strain a linear combination of tSZ and kSZ, with the 
resuh Dlf^^o + 0.5 D^^^^ = 4.5 ± 1.0 /xK^. For the base- 
line model with the kSZ power spectrum fixed to that 
expected for homogeneous reionization, we constrain the 
amplitude of the tSZ power at ^ = 3000 to be D^^^q = 
3.5 ± 1.0 ^K^. We test the dependence of the measured 
power on the specific tSZ model power spectrum and 
find that with the present data, the result is insensi- 
tive to the model shape for several recently published 
models. We also find that the measured tSZ power is 
insensitive to the modeling assumptions for the DSFG 
and radio foregrounds. We derive upper limits on sec- 
ondary anisotropy power when the kSZ amplitude is al- 
lowed to vary. At 95% CL, we find < 5.3 /iK^ and 



The corresponding limit on the total 
0.5D^^i?n < 6.2 at the 95% CL. 



^3000 

SZ power is D^fl, -r u.^ ^3000 

We compare the measured amplitude of the tSZ power 
spectrum with the predictions of recent models and sim- 



ulations. LIO found less tSZ power than the lSehgal et al., 
( 2010[ SIO) model predicts; the results in this paper in- 
crease the significance of that claim. We find that, for 
the cosmological parameters preferred by the primary 
CMB power spectrum, the simulations of SIO overpredict 
the tSZ power spectrum amplitude, with an observed-to- 
predicted ratio of Atsz = 0.33 ± 0.15. The more recent 



models of Shaw et al. 



(2010) and Trac et al. (20101 are 



gives Atsz = 0.59 ± 0.25. All three of these models as- 
sume that intra-cluster gas resides in hydrostatic equilib- 
rium in the host dark matter gravitational potential, and 
include prescriptions to account for st ar form ation and 
AGN feedback. The iTrac et all (|20T0| and Shaw et al. 



less discrepant, but still overestimate t he tSZ signal. For 



example, the 'nonthermal20' model of Trac et al. (20101 



(20101 models also include a significant non-thermal con 
tribution to the total gas pressure, which results in a 
lower tSZ signal. For the cosmological hydrodynamical 
si mulations including star formation and AGN feedback 
of lBattaglia et al.| p.010 ), we find Atsz = 0.45 ± 0.20. 
If instead we adopt the tSZ template derived from their 
simulations that do not include star formation and AGN 
feedback we find Atsz = 0.27 ± 0.10 and thus a larger 
discrepancy between the measured and predicted signal. 

We allow for a correlation between the tSZ and DSFG 
components in the model to investigate whether such a 
correlation could be responsible for the measurement of 
lower than expected tSZ power. With the kSZ ampli- 
tude fixed, marginalizing over the tSZ-DSFG correlation 
coefficient has no significant effect on the resulting tSZ 
power. When the kSZ power and tSZ-DSFG correlation 
are both free to vary, the fits actually prefer a larger kSZ 
power and slightly reduced tSZ power. We also note that 
if tSZ-DSFG correlation did significantly affect the mea- 
sured tSZ power, one would expect a difference between 
the result of the multi-frequency fits that assume no tSZ- 
DSFG correlation and fits to the differenced bandpowers, 
which should be largely free of DSFG signal. Instead, 
both analyses find the same tSZ power. We therefore 
conclude that correlation between the tSZ and DSFG 
components cannot be used to reconcile the measured 
tSZ power with model predictions, at least when oper- 
ating under the reasonable assumption that the Poisson 
and clustered DSFG spectral indices are similar. 

We also investigate combining the measured tSZ am- 
plitude with CMB data to obtain joint constraints on 
(78 ■ We find that, if we assume that the template predic- 
tions are perfectly accurate, the statistical uncertainty 
is reduced by as much as a factor of two compared to 
the primary CMB-only constraints. The resulting val- 
ues of as range fro m 0.771 zt 0.0 13 for the SIO model to 
0.802 ±0.013 for the ITrac et al.| ( r2010^ predictions. If we 
assume that the predicted amplitude is uncertain at 50% 
in order to account for the range of model predictions, 
the tSZ measurement adds little additional information, 
and the constraint degrades to as = 0.812 ± 0.022. Cur- 
rent modeling uncertainties prevent the use of the tSZ 
power spectrum as a tool for precision cosmology, how- 
ever, precise measurements of the tSZ power spectrum 
can provide new information about cluster physics in the 
poorly studied regimes of low mass and high redshift. 

In addition to the constraints on SZ power and erg, 
these data allow us to place constraints on foreground 
components. The measured bandpowers are consistent 
with a simple foreground model consisting of clustered 
DSFGs with a single spectral index and an unclustered 
distribution of radio sources (with a strong prior from 
radio source count models). Using this model, we find 
the spectral index of the DSFGs is a = 3.58 ± 0.09. The 
Poisson amplitude is -D3000 = 7.4 ± 0.6 /iK^ and the clus- 
tered amplitude is D^qqq = 6.1 ± 0.8 ^K^. We explore 
a number of variations on this simple DSFG model and, 
in general, find these variations have little impact on the 
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SZ constraints. By comparing fits using alternate models 
of the clustered DSFG power spectrum, we find the data 
show a slight preference for the power-law model over the 
linear-theory model used in HIO. 

We find a Poisson DSFG spectral index of = 
3.58 ± 0.09. Although the Poisson spectral index is well- 
constrained and fairly robust to changes in the model, 
the amplitude of the Poisson power changes by 4 cr as we 
change the shape of the clustered template. The Poisson 
power estimates are therefore dominated by this model- 
ing uncertainty, and any analysis of these results needs 
to take that uncertainty into account. The constraint on 
the clustering spectral index has improved from 3.8 ± 1.3 
reported in HIO to = 3.47 ± 0.34. This spectral index 
is sensitive to the shape of the clustered DSFG template: 
with the power-law template, we find ac = 3.79±0.37. In 
either case, the data are consistent with identical indices 
for both the Poisson and clustering DSFG contributions 
to the power spectrum. 

In 2009, the SPT receiver was refurbished to add sen- 
sitivity at 95 GHz. Since then, we have mapped an addi- 
tional 1300 deg^ to depths of 42, 18, and 85 /iK-arcmin at 
95, 150, and 220 GHz respectively. These data will break 
the tSZ/kSZ degeneracy present in the current data and 
facilitate the separation of the SZ components. Improved 
measurements of the tSZ power spectrum will provide 
new insight to the cluster astrophysics relevant in the 
poorly understood regimes of low mass and high redshift. 
We also expect these data to provide new constraints on 
the kSZ power spectrum, which have the potential to 
illuminate the epoch of reionization. 
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